多条人工地震测深剖面资料联合反演首都圈三维地壳结构

为了对多条二维剖面资料进行统一分析和解释,本文考虑了各剖面交点处的约束条件,通过建立统一的偏导矩阵,将各剖面资料联立求解.利用该方法对中国地震局地球物理勘探中心在首都圈内6条人工地震宽角反射/折射剖面资料进行了重新处理和解释,得到该地区的地壳三维速度结构和构造特征.在此基础上,结合其他地球物理研究成果,确定了壳内存在的4处低速层（体）的空间范围及断裂的走向.     

基于螺旋线上谱因式分解的地震波场隐式辛算法

均匀介质、复杂各向同性介质和各向异性介质中的地震波传播过程,可用统一形式的标量声波方程描述.考虑到在无损耗条件下,地震波方程描述了地震波场这一个无穷维的哈密顿体系随时间的演化过程,该过程为一个单参数连续辛变换,因而可以在其哈密顿形式表述下导出其辛格式.与显式辛算法相比,隐式辛格式对应的隐式辛几何算法具有无条件稳定的特点,可以允许较大的计算步长.但是由于隐式算法不可避免地面临高阶矩阵的求逆,其每一步的计算速度较慢.为实现矩阵快速求逆,文中采用了螺旋边界条件下谱因式分解的方法.在螺旋边界条件下,需要求逆的矩阵化为带状矩阵,而且其各列非零元素的位置和大小具有非常好的相似性,因而可以采用谱因式分解的方法实现快速LU分解.文中采用二阶精度的隐式蛙跳辛格式和谱因式分解方法,计算了常速度、层状介质和Marmousi模型中的波场.计算表明,隐式辛算法不失为波场计算的一种好方法.     

苏北大陆科学钻探靶区重磁异常反演解释

对江苏东海县毛北地区重、磁观测异常进行了处理和特征分析,并根据科学钻探主要目标体榴辉岩与其围岩有达0.9×103kg/m3的明显密度差,进行了等效密度反演成像与综合地震约束反演计算.给出了该地区重力异常模型,可作为综合地质与地球物理研究的一个基础.     

水与流体饱和孔隙介质界面上非镜面反射声场的实验研究

用光学法和发射-接收法研究了水与流体饱和孔隙介质界面上声波的非镜面反射声场．测量了入射角θ_i为16°－80°范围内的反射系数R，记录了θ_i在Rayleigh波临界角θ_c附近±4°范围内两束反射波的反射角、波束位移△、波形和频谱．实验发现：(1)在θ_c附近±4°范围内，观察到θ_i两束反射波A和B，其中波束B沿界面有波束位移△，实验测出△随θ_i增大而减小，并非在θ_i=θ_c时最大．(2)反射角等于θ_i的波束的频谱很接近于入射波的频谱，反射角等于θ_c的波束的频谱相对于入射波有频谱的负偏移．(3)波束A的相位与发射换能器直达波的相位一致，波束B的相位与波束A相反     

1998年5月磁暴磁层电流体系的地磁效应分析

低纬度地区地磁场的短时变化主要由以下电流体系产生:电离层发电机电流(IDC)、对称环电流(SRC)以及由部分环电流和Ⅱ区场向电流及其电离层回路组成的内磁层三维电流体系(PRFI).此外,由Ⅰ区场向电流及其电离层回路组成的电流体系(FACI)所产生的低纬地磁场也是不可忽略的.本文针对1998年5月1-6日的大磁暴,首先利用多个同子午线台站对的数据分离并消去由IDC电流产生的Sq场.然后,通过线性建模分离其他电流体系产生的磁场成分.结果表明:(1)发生在5月1-6日的磁暴可以分为两个过程,PRFI和FACI电流体系在1-3日不明显,在4-5日伴随着亚暴强烈发生.(2)SRC的变化情况在第一阶段同Dst指数相似,在第二阶段明显滞后于Dst指数.(3)在5月4-5日,PRFI在SRC之前增强,随着PRFI和FACI的恢复,SRC开始增强.这一结果为我们了解环电流和场向电流的形成以及它们的关系提供线索.     

辽宁青城子铅锌银金矿田控矿构造与找矿方向

辽宁青城子铅锌银金矿田为一集铅锌银金多个大型矿床于一区的重要成矿集中区。青城子中深构造层次的片麻岩穹隆（或变质核杂岩）构成了该区构造的主体。矿区中由片麻岩穹隆所构成的构造组合控制了矿体的定位、形态以及矿床的规模。青城子片麻岩隆起西南侧的卷曲褶皱翼内北东、北西走向的陡倾斜断裂发育，矿床主要为沿断裂交代充填的脉状铅矿，矿床规模较小；片麻岩隆起滑脱伸展翼多层顺层滑脱断裂带，为后期成矿提供了较大的容矿空间，有利于形成规模较大的矿床。     

酸性土壤中20种元素有效态浸提体系浅析

通过对酸性土壤中Cu、Pb、Zn、Co、Ni、Fe、Mn、Cd、K、Na、Ca、Mg、Si、P、S、Mo、F、Se、As、Hg20种元素有效态的浸提与测定，拟定了一个操作条件宽松、质量稳定，能满足区域调查需要的浸提体系和分析方案。该分析配套方案对20种元素的检出限在（0.0001～4.9）×10     

地震CT在寻找隐伏铜镍矿中的应用

用透射地震ＣＴ方法在金川龙首矿进行了隐伏铜镍矿体的探测试验．试验是在地下５００ｍ的１２８０平面上进行的，覆盖面积３００ｍ×３００ｍ．使用ＥＳ－２４０１Ｘ型４８道高分辨数字地震仪．铜镍矿体为高波速值的超基性侵入岩，采用纵波速度作为层析物理量是敏感的．研究区划分成８１×６４个单元，单元体大小为４ｍ×４ｍ．有效地震射线７６１９条，使用联合迭代重建法（ＳＩＲＴ）求解慢度走时方程．速度图像及地质解释表明测区中西部有一个高波速异常体，推断与含铜镍矿体或花岗体有关；此外还发现一条近东西向的低波速带，推断为断裂．     

非单一电离层扰动场的台阵探测分析

真实的ＴＩＤ一类大尺度电离层扰动是非平稳和色散的，而且是统计非单一的．本文首先给出了描述这种真实扰动场的数学模型，进而导出了联系台阵探测可测量与扰动场传播参量的观测方程，在此基础上提出了求解观测方程的时频分析方法．这一台阵探测数据分析的新方法，既能分离出非单一电离层扰动场中各种准单一扰动成分，又能获得准单一扰动成分传播参量随时间的演变和随频率的色散．分析实例表明，本文方法分辨率好，精度高，结果可靠，是适用于ＴＩＤ一类大尺度电离层扰动台阵探测数据分析的有效方法．     

射线追踪的微变网格方法

本文给出一种适用于叠前数据速度分析和叠前深度偏移的快速射线追踪方法──微变网格法；该方法不仅精度高，而且计算速度极快，可适用于当前计算条件下的地球介质速度反演和叠前深度偏移成像的要求．