垂向非均匀介质中转换波转换点计算

转换点位置的计算是转换波资料处理中的一个关键问题. 本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法. 研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法. 结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置.     

有限元方法在地震孕育过程与前兆机理研究中的应用进展

有限元方法的基本思想是将连续的求解区域离散为一组有限个、 且按一定方式相互联结在一起的单元的组合体。 有限元方法是目前应用最为广泛的数值模拟方法, 由于能够将复杂介质的力学本构关系及边界条件问题转化成为常规问题的计算程序, 所以越来越多的地震学者将其应用于地震成因与前兆机理研究中。 该文着重介绍了有限元数值分析方法及其模拟在地震孕育过程与前兆机理研究中的应用进展, 并提出开展小尺度规模如岩石中矿物颗粒之间的力学行为的有限元模拟研究等几点建议。     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.     

再论纬向正常密度假说与内波假说的比较

本文对“纬向正常密度假说”和“内波假说”再次进行了讨论,进一步比较了二者的联系和区别.根据岩石圈板块概念和地幔对流假说的基本观点、以及“纬向正常密度假说”和“板块运动重力机制”的基本观点,提出了“软流圈板块”概念.     

波形差异反演方法及其在气藏描述中的应用

在油气藏开发中后期,如何利用的已有的勘探成果和丰富的油气藏知识,加深对勘探区域油气开发潜力的认识,为油气田的稳产高产服务,是摆在油气田中后期开发人员面前的难题.油气田开发的研究人员,利用现有的资料,开发出了很多方法,结合地震测井等资料为油气田的稳产高产做出了巨大的贡献.本文在分析了地震波在地质体中传播过程中能量变化,相应的产生了振幅变化的理论,利用三维地震数据,通过对三维地震数据解释层位对应的波形差异的分析,结合工区内的生产井和开发井的对应地层中的油气信息,进行反演得到地层的波形差异数据体,利用相应的绘图软件绘制二维平面图形,结合工区内的沉积相等地质信息对其数据体进行油藏解释,最终得到波形差异反演数据体,利用波形差异反演数据体进行油气藏的识别、描述和预测,为油气田的中后期开发服务.     

完全匹配层吸收边界在孔隙介质弹性波模拟中的应用

模拟弹性波在孔隙介质中传播,需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上,首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先,本文通过均匀孔隙模型的数值解与解析解的对比,验证所提出的数值方法的正确性;然后,本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能,将完全匹配层与廖氏和阻尼吸收边界进行了对比,研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明,在孔隙介质中,完全匹配层作为吸收边界能十分有效地吸收衰减外行波,无论对体波还是面波,是一种高效边界吸收算法.     

北太平洋对人为二氧化碳吸收的数值模拟

采用一个开边界海盆尺度环流模式研究人为CO2在北太平洋的吸收和分布,并与闭边界模式的模拟结果进行了比较.模拟结果表明,西北太平洋和赤道东太平洋是两个重要的人为CO2汇.使用较大的等密度面扩散系数使得西北太平洋人为CO2通量增大,而赤道海区通量减小（RUN2）.与闭边界模式相比,开边界模式中该两个区域的人为CO2通量都增加了.1800～1997年间,北太平洋共吸收人为CO2 2375GtC(1Gt=1×1015 g, RUN1).副极地海区是人为CO2的一个重要输出区,能输出人为CO2吸收量的38％～54％,而20°N～30°N海区是人为CO2的一个重要贮存区,占整个北太平洋的24％.开边界对于10°N以南太平洋吸收和贮存人为CO2有很大影响.与基于观测资料的估计相比,虽然模式低估了西北太平洋的人为CO2的穿透,高估了东北太平洋的人为CO2的穿透,但总的说来,模式比较好地估计了人为CO2在北太平洋的贮存.     

使用纵波和横波走时层析成像研究菲律宾海板块西边缘带和南海地区的深部结构

使用纵波震相及同台同源纵波和横波震相的走时资料及层析成像反演方法,分别给出了菲律宾海板块(PHP)和南海地区的纵波速度V_pf及同台同源的纵波和横波速度V_p和V_s结构. 结果表明,（1） PHP与欧亚板块(EUAP)的俯冲接触关系随地段而异,在琉球海沟,PHP向EUAP之下俯冲深达400 km;在台湾岛,EUAP先近陡直地俯冲到深度400 km,然后折向PHP之下达到660 km左右;在马尼拉海沟北段,俯冲板片几乎垂直达到660 km附近;在菲律宾海沟,PHP向EUAP之下的俯冲深度超过660 km．（2）南海地区之下是一个深达300～400 km的宽阔低速区,并且大致在莺琼海盆与700 km深处另一低速区曲折相通;在该宽阔低速区内部,有两个小而明显的低速区分别在海口火山和珠江口盆地下方．（3）对V_p和V_s及据其算出的容声速度V_b作分析发现,V_s和V_b的平均扰动量对深度的变化在一些深度范围内是反向的;年龄较大的太平洋板块俯冲体的V_s相对扰动量的平均值大于V_b的,而在较年轻的PHP俯冲体中则正好相反．     

中国地震台网震级的对比

采用线性回归和正交回归方法,利用中国地震台网1983-2004年的观测资料,对中国地震局地球物理研究所测定的地方性震级ML、面波震级MS与MS7、长周期体波震级mB、短周期体波震级mb进行对比,给出了它们之间的经验关系式.研究结果表明:①由于不同的震级标度反映了地震波在不同周期范围内辐射地震波能量的大小,因此对于不同大小的地震,使用不同的震级标度更能客观地描述地震的大小.当震中距小于1 000 km时,用地方性震级ML可以较好地测定近震的震级.当地震的震级M＜4.5时,各种震级标度之间相差不大.当4.5＜M＜6.0时,mB＞MS,MS标度低估了较小地震的震级,因此用mB可以较好地测定较小地震的震级.当M＞6.0时,MS＞mB＞mb,mB与mb标度均低估了较大地震的震级,用MS可以较好地测定出较大地震(6.0＜M＜8.5)的震级.当M＞8.5时,MS出现饱和现象,不能正确地反映大地震的大小;②在我国境内,当震中距＜1 000 km时,ML与区域面波震级MS基本一致,在实际应用中无需对它们进行震级的换算;③虽然MS与MS7同为面波震级,但由于所使用的仪器和计算公式不同,MS比MS7系统地偏高0.2～0.3级;④对于长周期体波震级mB和短周期体波震级mb,虽然使用的计算公式相同,但由于使用的地震波周期不同,对于mB=4.0左右的地震,mB与mb几乎相等,而对于mB≥4.5的地震,则mB＞mb.     

共炮检距道集波动方程保幅叠前深度偏移方法

本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息.