一个中纬电离层E层理论模式

本文建立了一个中纬电离层E层理论模式. 该模式从NO+,O２＋,O+和N２＋这四种主要离子的连续性方程出发,通过数值模拟得到中纬电离层E层电子和各种离子密度随时间和高度的变化情况. 计算结果能较好地反映出E层电子密度峰值（NmE）或E层临界频率（foE）的日变化、季节变化以及随太阳活动的变化趋势. 将模式的计算结果与武汉地区测高仪的观测数据进行比较,结果证明模式能够较为客观地反映中纬电离层E层的实际形态. 针对以往电离层E层理论模式存在的主要问题,本文还进一步讨论了几种重要因素,包括二次离化源,λ<150?谱段的辐射通量,吸收截面以及NO分布对于模式计算结果的影响.     

利用Cluster Ⅱ数据分析等离子体片区磁结构运动特性

利用Cluster Ⅱ 卫星上搭载的磁通门磁强计的观测数据,可以计算磁结构的运动速度.本文在GSM（Geocentric Solar Magnetospheric System,地心太阳磁层系）坐标系下处理了磁暴期间（2002年8月1日到2日）和磁场宁静期（2002年8月6日到7日）的等离子体片区的观测数据.通过分析磁结构速度的统计特性和频谱特性,得到以下结论：（1）磁结构运动速度呈现出不规则的摆动结构.（2）Vz的方差与Vy的方差近似相等,而远大于Vx的方差.这可能是由于“风袋效应”.（3）在黎明侧的等离子体片区, Vx和Vz的平均值约为0,但Vy的平均值的量级为-10 km/s.这可能是在等离子体片区Pc5压缩波造成的.（4）Vx, Vy, Vz的频谱在0033 min-1和0035 min-1之间都有一个峰值.     

基于二阶矩闭合模式平衡解的海浪生湍流混合系数解析估计

从次小尺度海洋运动方程和运动近各向同性假定出发,本文给出高确定性的海洋湍流二阶矩闭合方程和特征量变化方程.这组方程描述波动、涡旋和环流等较大尺度运动对海洋湍流的剪切不稳定性生成机制.同时,基于海浪破碎是湍流海面生成主要机制的认识,本文按破碎海浪统计理论给出湍流特征量变化方程的海面边界条件,包括由破碎能量损耗和卷入深度分别确定的海面动能通量和动能耗散率.根据上层海洋湍流特征量方程的平衡解和关于动能耗散率的观测资料,本文得到海浪生湍流混合系数的解析估计.所得到的海浪生湍流动能耗散率与观测资料进行了定性与定量对比,混合系数与先前Prandtl混合长度理论估计的一致性检验.     

对橄榄石相变形核率的估算

地幔主要矿物橄榄石［(Mg089Fe011)2SiO4］的相变与深源地震的解释以及俯冲带和周围地幔的相互作用相关. 其中，形核率和长大率是刻画其相变动力学的两个重要参数. 因为实验技术的原因，目前人们还没有基于实验数据给出橄榄石的形核率. 本文首先通过对挤碰物理图像的分析，对边界形核情形的相变动力学理论中表征挤碰程度的无量纲参数进行了修正. 在此基础上，根据已有的实验资料，首次对橄榄石相变的形核率进行了估计. 通过对形核率数据的拟合，可以得到形核率参数K0和γ1/3. 得到两个形核率参数的变化范围分别为K0=55×1021～87×1027s-1·m-2·K-1，γ1/3=0～020 J·m-2，最佳拟合值为：K0=69×1024s-1·m-2·K-1，γ1/3=016 J·m-2.     

青藏高原东北缘地壳三维速度结构

本文用1980—2000年M≥1.5的2 032个天然地震事件的38 052个〖AKP-〗、〖AKS-〗、Pm、Sm、Pn和Sn震相到时及人工地震测深给出的Moho面形态资料,利用地震层析技术反演了32°~40°N, 100°~108°E区域内地壳地震波速度结构.从层析成像图象中可以得到,本区的地壳可分成4个层位.第1层（埋深约在0~3 km）为沉积层, 速度梯度约为0.2 s-1;第2层（埋深约在3~17 km）为上地壳, 其顶部速度梯度约为0.1 s-1, 下部速度横向变化较大且存在低速块体;第3层（埋深约在17~36 km）为中地壳, 速度梯度约为0.03 s-1;第4层（埋深约在36 km—Moho）为下地壳, 是一个契形层,总的趋势是西厚东薄,青藏高原较厚逐渐向鄂尔多斯地块和扬子准地台方向变薄,各处的地震波速度梯度不尽相同.     

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)的重大影响.     

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)的重大影响.     

岩石的纵、横波速度与密度的规律研究

选择了济阳坳陷第三系地层四类不同岩性、237块典型岩心样品,在变温变压条件下进行超声波纵、横波速度及密度的测定,样品的孔隙度范围从2％到27％.根据变温变压测试数据,得出纵横波速度的压力、温度校正公式,将高温高压条件下测量的弹性参数校正到60 MPa条件,再进行密度、Vp和Vs数据分析.分析结果表明：Vp、Vs具有很好的线性相关性,得到新的密度与Vp关系式.同时将Vp、Vs关系融入到Gardner公式中,将Gardner函数推广为〖JP2〗ρ=CVApVBs的函数形式（称为广义Gardner公式）,拟合出砂岩的密度、Vp和Vs间的新关系式：ρ=1.575 8V0.287 9pV0.019 7s.〖JP〗文中讨论了该公式在横波速度预测、岩石物理模拟和AVO正反演等的作用,指出这些规律性认识可服务于油气勘探和开发.     

波浪破碎卷入气泡的泡径分布理论模型

借助相似性定理推导出波浪破碎形成的气泡总数与气体体积卷入率、湍流能谱密度以及表面张力之间的关系, 得出泡径谱的指数为-2.5+n_d, 其中n_d取决于湍流能谱在黏性耗散区的分布. 利用引入的两个破碎统计量, 进一步导出了破碎过程中的气体体积卷入率Q的表达式, 建立了泡径谱与波浪要素的联系. 在此基础上, 利用观测结果提出两点假设, 并据此将泡径谱N(a)推广为随泡径和深度的分布函数N(a, z).     

微粒对多晶冰流变行为的影响——(Ⅱ)滞弹性行为

用反相直接加载的方法研究了微粒对多晶冰滞弹性的影响.结果表明：在较高频率时(1 Hz 和 10-1 Hz)，滞弹性应变与应力峰值呈线性关系；在较低频率时(10-2 Hz)，其应力/应变偏离线性关系.微粒在高频时(1 Hz)通过阻碍晶界滑移对晶界弛豫产生重要影响，增加了模量和降低了内耗.但微粒对低频时的位错弛豫没有明显的影响.通过滞弹性实验计算了非线性弛豫模型的两个重要参数，计算表明K值约为007 Pa，α值约为054.     

基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼Mw75级地震参数－Ⅰ.均匀滑动反演

1997年11月8日西藏玛尼Mw75级地震是干涉雷达技术应用于地震观测以来的一次重要事件.Peltzer等认为该形变场体现了地壳的非线性弹性形变响应，挤压象限和拉张象限的杨氏模量比值为2.相同的情形并没有在其他地震同震形变场中发现，通常用线弹性理论就能够很好地解释观测数据.考虑到他们所用断层模型的简化程度和纯走滑约束等，本研究认为非线性弹性解释是牵强的.应用广泛使用的Okada线弹性位错模型，采用卫星观测得到的地表断层位置，去除倾滑为零的约束，基于非线性优化反演方法寻求拟合雷达观测的最佳断层几何参数和均匀滑动参数.结果表明，线弹性模型能够满意解释观测数据.断层在朝阳湖以东的段落最深达到165 km，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾     

2002年南海季风爆发前后西沙海区海-气通量交换及其变化

利用2002年4月24日至6月20日在西沙海区进行的第三次南海海-气通量观测试验资料,采用涡相关法和TOGA COARE25b版本通量计算方案,计算了西南季风爆发前后海洋-大气间的通量交换,讨论了辐射、动量、感热通量、潜热通量、海洋热量净收支的时间变化特征及其与气象要素变化的关系.结果表明：西南季风爆发前后,太阳短波辐射、海面净辐射、潜热通量和海洋热量净收支变化特别强烈;通量变化受不同环境要素的影响：感热通量与海-气温差呈正相关关系,与气温呈明显的负相关关系.潜热通量与风速、海-气温差及海面水温均有正相关关系,其中与风速的关系最密切.动量通量（τ）主要随风速变化,它与风速（V）的关系可以表示为τ=000185V2-000559V+001248.     

云闪放电对云内电荷和电位分布影响的数值模拟

我们应用改进的随机闪电参数化方案，对两次雷暴的起、放电过程进行了二维125 m和250 m分辨率云闪模拟试验，分析表明高分辨率模拟的云闪通道几何结构、延伸范围和最大垂直电场变化等特征与观测结果更为一致，并且揭示了云闪放电重新配置云内电荷分布和空间电荷中和过程的一种新的物理图像：(1) 云闪的直接物理效应是在已有的空间电荷区内沿着放电通道沉积异极性电荷、形成复杂的空间电荷分布，有效地导致云中电位和场强绝对值及静电场能量剧烈下降，使放电终止.其中：放电后在闪电通道经过主要区域，电位降到±30 MV之间，垂直电场强度也降到±20 kV·m-1之间，一次正或负云闪估计消耗掉107～1010 J静电能；(2) 云中电荷中和不是正、负空间电荷简单地直接相互抵消的瞬态过程，而是广泛分布的云中空间电荷与通道沉积的异极性电荷通过湍流交换、平流输送、重力沉降以及起电等多种因素逐渐混合并部分被中和的一个后续慢过程，其弛豫时间典型值在14～44 s之间，在此期间通道感应电荷总量下降到50%以下.并残余部分电荷参与后续放电前高空间电荷密度和强电场的重建过程.     

一次磁层亚暴期间磁场偶极化的尾向传播

利用TC1、Cluster和Polar结合极光和同步高度及地磁的观测，研究了2004年9月14日1730～1930 UT时间段的亚暴偶极化过程.此前行星际磁场持续南向几个小时.亚暴初发(Onset)开始于1823 UT.2 min之后，同步高度的LANL 02A在子夜附近观测到了明显的能量电子增强（Injection）事件，而TC1在1827UT左右在磁尾(-10,-2, 0)RE (GSE)观测到了磁场BX的突然下降，伴随着等离子体压强和温度的突然增加及磁场的强烈扰动.在(-16, 1, 3)RE (GSE) 的Cluster上相同的仪器观测到相同的现象，只是比TC1观测到的晚大约23 min，在1850 UT左右.虽然Polar在更靠近地球的较高纬度(-75, 35, -40)RE (GSE)附近，也在1855 UT左右观测到了这种磁场偶极化现象.以上的观测时序表明TC1、Cluster观测到的磁场偶极化比亚暴偶极化初始发生分别晚4 min和27 min.说明偶极化由近磁尾向中磁尾传播.详细计算表明偶极化源区的位置大约在X=-77RE～-86RE，而传播速度大约为70 km·s-1.在这个事件中亚暴的物理图像可能是中磁尾的近地重联产生的地向高速流到达近磁尾，为近磁尾的亚暴触发创造了条件；亚暴在近磁尾触发之后，磁场偶极化峰面向中磁尾传播.     

阿尔奇(Archie)公式：提出背景与早期争论

根据关于岩石电学的早期资料，介绍了阿尔奇公式的提出背景和当时身处学术界和测井界的研究者和工程师们围绕阿尔奇公式所展开的一系列研究和论战.回顾与分析证明：围绕阿尔奇公式所进行的早期论战在于地层因子F及指数m和n.具体的论战焦点有四个：（1）地层因子的表达式F=1/m是否合适；（2）m和n是否具有物理内涵；（3）m和n的值到底是常数还是变数；（4）m和n与岩石的内部结构是否有关.在经过近20年的研究和论争后，即在1960年左右发现指数m和n与岩石的结构、构造以及孔隙流体的饱和度和流体替换的历史有关.此外，还发现岩石的湿润性和超毛细管孔隙度对岩石的电阻率也有很大的影响.虽然阿尔奇公式影响深远而且在岩石物理学和石油测井的发展中具有划时代的意义和里程碑式的作用，但是没有考虑泥质、淡水、低孔隙、非均匀几何参数分布（孔隙度、曲折度）、非均匀饱和度分布、各向异性以及参数a和b对岩石电阻率的贡献和隐藏在这些参数和因素后面的物理机制.因此，在应用阿尔奇公式时必须注意有关的应用条件.此外，根据完全水饱和公式ρ0=FρW=aρW/m和部分水饱和公式ρt=Iρ0=bρ0S-nW来研究F、I、a、b以及m和n与岩石的内部结构、孔隙分布以及孔隙流体和岩石颗粒之间的物理及化学相互作用的关系具有反问题的特点.所以，围绕阿尔奇公式的研究和争论将有可能会伴随着时代的进步和科学技术的发展一直进行下去.