马蹄形隧洞考虑支护滞后过程的应力分析

基于平面弹性复变函数中的保角变换方法,推导出带有衬砌的非圆形隧洞在原始地应力作用下的应力解析解。根据衬砌内边界的应力边界条件及围岩衬砌接触面上的应力和位移连续条件,获得求解围岩和衬砌解析函数的基本方程,计算了围岩和衬砌中的应力和位移。在求解过程中,考虑了支护滞后于开挖的力学过程,并认为围岩和衬砌之间紧密接触,不会相互分开和相对滑动。以马蹄形隧洞为例,获得了围岩开挖边界和衬砌内外边界的切向应力及围岩与衬砌接触面上的接触应力分布规律,并与ANSYS数值方法结果对比,算例表明两种方法的计算结果吻合很好。     

基于加载动力效应对孔压滞后现象的一个解释

在底部可测孔压的固结仪内进行单面排水固结试验时,观测到试样底部的孔压可出现滞后现象。为分析其原因,考虑加载引起的动力效应,建立以竖向应变为求解对象的一维动力固结方程,并用有限差分法进行求解。引入新指标CD来综合考虑试样压缩模量、渗透系数、密度以及土样厚度等4个参数对动力效应的影响,探讨了该指标对固结过程中孔隙水压力的影响。计算结果表明,瞬时加载引起的动力效应是引起试样内孔压滞后的一个原因,并且对于压缩模量或渗透系数较大的试样,更易出现孔压滞后现象,同时动力效应对固结度的影响则主要体现在固结的前期。最后讨论了固结分析时可以忽略动力效应的范围。     

磁暴期间热层大气密度变化

基于CHAMP卫星资料,分析了2002—2008年267个磁暴期间400km高度大气密度变化对季节、地方时与区域的依赖以及时延的统计学特征,得到暴时大气密度变化的一些新特点,主要结论如下:1)两半球大气密度绝对变化(δρa)结果在不同强度磁暴、不同地方时不同.受较强的焦耳加热和背景中性风共同作用,在北半球夏季,中等磁暴过程中夜侧和大磁暴中,夏半球的δρa强于冬半球;由于夏季半球盛行风环流造成的扰动传播速度快,北半球夏季日侧30°附近大气,北(夏)半球到达峰值的时间早于南(冬)半球.而可能受半球不对称背景磁场强度所导致的热层能量输送率影响,北半球夏季强磁暴和中磁暴个例的日侧,南半球δρa强于北半球;春秋季个例中日侧30°附近大气,北半球先于南半球1~2h达到峰值.2)受叠加在背景环流上的暴时经向环流影响,春秋季暴时赤道大气密度达到峰值的时间最短,日/夜侧大气分别在Dstmin后1h和2h达到峰值.至点附近夜侧赤道大气达到峰值时间一致,为Dstmin后3h;不同季节日侧结果不同,在北半球冬季时赤道地区经过更长的时间达到峰值.3)日侧赤道峰值时间距离高纬度峰值时间不受季节影响,为3h左右.在春秋季和北半球冬季夜侧,赤道大气密度先于高纬度达到峰值,且不同纬度大气密度的峰值几乎无差别,表明此时低纬度存在其他加热源起着重要作用.     

基于BDS与GPS精密单点定位的天顶对流层延迟估计比较

在BDS与GPS现有星座条件下,针对若干IGS和MGEX跟踪站的实测数据,利用CODE事后GPS产品与WHIGG计算的BDS精密轨道和钟差,对GPS单系统、BDS单系统及两者组合系统进行精密单点定位（PPP）处理,估计出相应的天顶对流层总延迟量,并进行分析比较。实验表明,与IGS提供的对流层产品相比,利用GPS单系统处理,能较准确地反映出天顶对流层延迟量,其精度为mm级;BDS单系统结果较GPS单系统略差,其精度优于2cm;GPS与BDS组合系统的结果与GPS单系统结果相近。     

电离层三维层析成像的自适应联合迭代重构算法

在电离层层析成像过程中,联合迭代重构算法是一种常用的反演算法.然而,该算法迭代收敛较慢,反演结果精度不高.为此,本文发展了一种自适应的联合迭代重构算法,该算法利用上一轮的电离层电子密度反演结果,自适应地调整松弛因子和加权参数.通过模拟数据和实测数据对该算法的反演结果进行了验证,并将得到的反演结果与电离层测高仪数据进行了比较,结果表明,该算法能够有效地反演电离层电子密度,且反演结果精度优于常用的联合迭代重构算法.     

GPS共视中电离层时延计算方法比较研究

为了更好地计算GPS CV（共视）时间传递中的电离层时延值（它是影响CPS CV比对结果精度的主要因素之一）,介绍了当前3种电离层时延的计算方法,并以NICT（National Institute of Information and Communications Technology）单站GPS比对数据及NICT与NTSC（National Time Service Center）的GPS共视比对数据为例,分析比较了不同的电离层时延计算方法对GPS时间比对结果精度的影响。计算结果表明：利用双频实测电离层时延和利用ICS（International GPS Service）提供的TEC（total electton content）map计算的电离层时延对GPS CV比对结果修正后的精度,比利用电离层改正模型的时延对比对结果修正后的精度分别提高30％～40％和20％～30％。     

高原重要机场航班延误的气象因素分析

通过对四川航空公司2005～2006年高原特殊机场航班延误的气象资料分析,得出影响高原机场最大的气象要素为低云、低能见度以及顺风超标等。并针对航班延误比较多的高原机场进行详细的数据统计以及气象特征分析,为排班工作提供参考。     

数值模拟静态应力扰动下的断层失稳:结果分析兼与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模型则不能完整地给出受到静态剪应力扰动后断层失稳时间提前或推后的估计值.     

对流层延迟模型映射函数研究

在介绍经验映射函数NMF和动态映射函数VMF1、GMF模型的基础上,利用ECMWF提供的格网数据,具体研究和分析了这3种映射函数1 a中的时间变化特征及其随高程、纬度、高度角变化的特点,并比较了3种映射函数随高度角的变化对斜路径延迟估计的影响。     

同时估计电离层延迟的单频精密单点定位方法

提出一种基于单频码和相位观测量的单频精密单点定位方法,将每个观测量的电离层延迟量与接收机钟差、对流层天顶延迟、接收机位置、相位模糊度一起作为未知参数。采用约化参数的平方根信息滤波与平滑算法进行参数解算。该方法适用于实时定位和事后处理,且不需要外部的电离层模型。采用全球分布的32个IGS监测站16 d实测数据进行静态解算试验,结果表明E、N、U方向的RMS分别为0.023 m、0.018 m、0.059 m;基于一组机载GPS数据进行动态解算试验,得到E、N、U方向的RMS(与载波相位动态相对定位结果比较)分别为0.168 m、0.151 m、0.172 m。