聂拉木地区高喜马拉雅岩石磁组构及其构造含义

利用岩石磁化率量值椭球体与岩石构造应变椭球体的共轴性,在缺乏岩性标志层的高喜马拉雅结晶基底（聂拉木地区）进行岩石磁组构研究,进而研究喜马拉雅造山带的变形期次和过程.磁组构特征表明：磁化率各向异性度P在樟木镇北1km处的片岩最高为172,在康山桥南1km处最低,为113;百分率各向异性度H在1182%～4507%之间;磁性线理L在102～109之间,磁性面理F在106～160之间,磁性面理F比磁性线理L发育;磁化率椭球体的形状因子T在019～076之间,磁化率椭球的扁率E在104～149之间,磁化率椭球体的形状为一压扁椭球;樟木镇北1km到肉切村,平均最小磁化率主轴方向D3为S－N（除聂拉木北1km的MA9外）,大多数样品最小磁化率主轴方向倾角I3＞51°,最高达726°,局部发育有拉长形磁化率量值椭球体.本文研究结果表明,聂拉木地区早期经历强烈的韧性变形,推测可能是一条巨大的右旋逆冲韧性变形带,晚期构造掀斜.     

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

速度差对井间走时层析成像质量的影响及其改进方法

在井间地震初至走时层析成像中,随着相邻地质体速度差的增大,使得射线分布不均匀,以及网格剖分不合适,导致层析成像结果不理想.物理和数值模型的井间走时层析成像表明：当速度差超过300%时,层析结果畸变较大;在300%～150%之间时,层析结果较好;低于150%时,层析结果好.在此基础上,提出了井间多尺度初至走时层析成像方法,即对同一模型采用多种网格剖分来同时进行层析成像,以获得研究区域的速度图像.数学和物理模型的井间多尺度走时层析结果表明：该方法很好地兼顾了层析成像的分辨率和精度,极大地改善了井间地震层析成像的质量.即使速度差超过30%,其多尺度的层析结果仍然较好.因此,这种方法具有实际应用的潜力.     

磁场重联中电子加速问题的物理分析和数值研究

观测表明耀斑中电子加速发生在软X射线耀斑环上方的磁重联区域.在重联电流片中被super Dreicer电场直接加速似乎是产生能量在10keV～10MeV之间高能电子的最直接的方式.本文的结果证明纵向电磁场可以有效地将电子“锁”在重联电流片上,使得横向电场得以直接加速电子.对于解释产生相对论性粒子的脉冲式耀斑,这可能是一个有效的机制.     

2004年3月18日高纬磁层顶多重通量管事件分析

2004-03-18 23:10～23:50 UT期间,“双星（Double Star）”探测一号卫星（TC 1）在向阳面磁层顶高纬晨侧由内向外穿越磁层顶,其时TC_1的GSM坐标为 (75RE, -55RE, -54RE), RE为地球半径.穿越过程中TC_1观测到了8个通量管和1个磁通量传输事件（FTEs）.在此期间Cluster星簇位于向阳面太阳风内,其GSM坐标为(180RE, -31RE, -62RE),其4颗卫星监测到行星际磁场（IMF）的BZ分量持续南向,BＹ有较大的负值.本文的研究表明：TC_1观测到的前7个通量管具有准周期重现性,周期大约是1～4 min,明显小于以前所观测到的FTEs的平均周期(8～11 min);所有的通量管都具有较强的核心场.本文分别使用最小方差分析法（MVA）和Grad_Shafranov反演方法（GSR）对通量管的轴向进行了分析和对比,发现所有的通量管主轴基本沿晨昏向,结果显示GSR方法在轴向分析上比MVA优越.本文使用GSR方法对通量管的磁场结构进行了分析,恢复出了通量管的磁场在卫星穿越面的结构图;此外,本文还对这次多重通量管事件进行了deHoffmann Teller(HT)分析,结果表明,所有通量管大致朝南极方向运动,均来源于向日面低纬区域.这说明它们可能起源于向日面低纬区,由该区的磁场分量重联产生.     

全球地幔垂直流动速度研究

用高分辨率地震体波速度成像以及相关的地球物理资料,计算地幔垂直流动形式及流动速度,得到全球地幔流垂直运动模式.从全球尺度来看,地幔流基本可划分为以下几个区域：欧亚大陆—澳大利亚、北美洲—南美洲为两个大规模下降流区域,西印度洋—非洲及大西洋、中南太平洋及东太平洋为两个大规模地幔上升流区域.地幔上升流起源于核幔边界,主要表现在地幔中部和上地幔下部.地幔垂直流动速度约每年1～4cm.地幔流动对地表板块运动、海洋中脊和中隆、俯冲带和碰撞带的分布起着控制作用.地幔上升流与地表现代热点有密切关系.从东亚尺度看,地幔流大体分为三个区域：东亚边缘裂谷系和西太平洋边缘海为上升流、西伯利亚地幔深度表现为物质下降流、青藏高原—缅甸—印度尼西亚特提斯俯冲带地幔下降流,这三个区域地幔流动与地表的西太平洋构造域、亚洲构造域和特提斯构造域相吻合.勾勒出南海地区构造特征：从上到下的大体结构是上部呈“工"字型、中间为圆柱型、底部呈盾形的地幔上升流.     

关于暴时电离层电流分布的南北半球不对称性

采用国际上广泛认可的高层大气和电离层经验模式提供的各种参数, 通过电离层电流连续方程, 计算出强磁暴条件下6月至日和12月至日内, 磁纬±72°和磁地方时00:00～24:00之间电离层电场、电流等的分布. 计算中考虑了地磁和地理坐标间的偏离; 除中性风场感生的发电机效应外, 还包含了磁层耦合(极盖区边界的晨昏电场和二区场向电流)的驱动外源. 结果表明, 6月至日时, 磁层扰动自极光区向中低纬的穿透情况在南、北半球内基本接近, 北半球内略强; 但12月至日时, 呈现明显的不对称性, 南半球的电流穿透远强于北半球, 而电场的穿透则是在北半球更强. 无论南北半球, 在中高纬地区, 午夜至黎明时段出现较强的东向电场分量, 其[WTHX]E×B[WTBZ]的向上漂移效应, 正是解释我们以往不少研究现象中所期盼的物理机制.     

地震雷达

利用绿色、环保的人工震源,主动向地下发射地震波,构建地震雷达,实现大范围的地下探测,是地球科学的一个前沿课题.最早,探测地球内部利用的是天然地震产生的地震波,这是因为天然地震释放的能量大,一次天然地震释放的能量相当于万吨级炸药.几千公里远处仍可以接收到信噪比大于1的地震波信号,但天然地震发生频度低和震源位置定位精度低限制了利用天然地震进行地下探测的精度;后来,利用人工源的地震勘探得到迅速的发展,探测精度明显提高,但受到人工源能量的限制,探测的空间尺度有限;今天,以小当量激发实现大尺度探测是发展地震雷达的一个关键问题,是将地震波理论和现代信息科学相结合的一个新的领域.发展人工震源的编码和接收信号的解码等理论和技术,可以从电磁波雷达在过去半个世纪所经过的道路得到许多借鉴和启发.可以预期,地震雷达的发展将会对观测地震学带来全面的影响.     

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

本文用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）为下地壳, 是一个契形层,总的趋势是西厚东薄,青藏高原较厚逐渐向鄂尔多斯地块和扬子准地台方向变薄,各处的地震波速度梯度不尽相同.     

河南省科研院所科技产业化融资问题研究

科技产业化是科研院所科技资源转化为经济效益、科研成果转化为现实生产力最直接的方式,而有效的融资是科技产业化的必备条件。河南省的科研院所是河南省科技创新的重要源泉。文章在提出该省科研院所科技产业化必要性的基础上,分析了科技产业化融资存在的问题、科研院所科技产业化所需的融资环境,并提出了解决融资问题的对策措施。