基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演

各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI(Message Passing Interface,消息传递接口)技术实现分频并行计算,测试结果显示并行接近线性加速比.     

地震预警现地PGV连续预测的最小二乘支持向量机模型

为提升现地仪器地震烈度预测的准确性与连续性,研究面向地震预警的PGV连续预测模型.以中国仪器地震烈度标准的计算参数:0.1~10 Hz带通滤波三分向矢量合成速度峰值PGV为预测目标,利用日本K-net与KiK-net台网P波触发后1~10 s强震数据,基于人工智能中的机器学习方法-最小二乘支持向量机,选取7种特征参数作为输入构建最小二乘支持向量机PGV预测模型LSSVM-PGV.结果表明,本文建立的LSSVM-PGV模型在训练数据集与测试数据集上的预测误差标准差变化趋于一致,具备泛化性能;P波触发后3 s预测PGV与实测PGV即可整体符合1∶1关系,随着时间窗的增长,PGV预测的误差标准差显著减小、并在P波触发后6 s趋向收敛,具备准确连续预测能力;对比同为P波触发后3 s的常用P d-PGV模型,LSSVM-PGV模型的PGV预测误差标准差明显减小,“小值高估”与“大值低估”现象明显改善,预测准确性得到提升.熊本地震序列的震例分析表明,对于6.5级以下地震,LSSVM-PGV模型最多在P波触发后3 s即可预测出与实测PGV整体符合1∶1关系的PGV;对于7.3级主震,由于其破裂过程的复杂性,P波触发后3 s的预测结果出现一定程度的低估,但随着时间窗增长至6 s时,预测PGV与实测PGV符合1∶1关系、并直到10 s整体趋势保持一致.本文构建的LSSVM-PGV模型可用于现地地震预警仪器地震烈度的预测.     

基于Levenberg-Marquardt法和最小二乘法的中国西部烈度衰减关系模型研究

收集1900年以来,中国西部地区488次地震的烈度等震线资料,采用长、短轴椭圆模型,应用Levenberg-Marquardt法和最小二乘法拟合了中国西南地区和西北地区地震烈度衰减关系,结果表明:两种拟合方法的回归模型结果具有良好的一致性,在西南地区与西北地区2个统计单元内地震烈度衰减关系有显著不同,应作为不同的分区对待.同时,将本文结果与其他研究者得出的中国西部地区地震烈度衰减关系结果进行了对比分析,并以近年来中国西部发生的6个真实地震的烈度圈数据为基准,对比了本文模型与其他区域模型的长、短轴及面积的计算结果,结果显示:本文烈度衰减关系计算结果的精准性要稍高于大区域的计算结果,但相对于小区域来说,本文结果比小区域的烈度衰减关系结果误差偏大.由于地震烈度衰减关系具有区域性特征,区域划分越细范围越小,结果精准性越高,但是计算大区域的烈度衰减关系,有利于第一时间掌握灾区大范围的灾情分布,特别是对于震级较大的地震,计算得到大尺度上的地震烈度衰减关系,可以迅速分析相邻省区受影响地区的烈度分布特征,无须利用另一个省的地震烈度衰减关系数据重新计算,故而本文两种方法得出的结果较好地反映了烈度衰减的区域特征,适合用于中国西部地区的工程地震研究和应用.     

衰减补偿的黏弹性最小二乘逆时偏移

地震勘探是发现和认识油气藏的重要方法,随着国内油气勘探领域、勘探对象的日趋复杂化,油气勘探迫切需求地震成像的精准度、分辨率和保真度的大幅度提升.相比于常规逆时偏移方法,最小二乘逆时偏移成像有较少的伪像,成像振幅均衡,成像分辨率高.目前常用弹性最小二乘逆时偏移方法来对地下介质进行成像.但由于地下介质通常表现出黏弹性性质,在含油气地层中表现尤为明显,若忽略地下介质的吸收衰减作用,地震波传播中的振幅衰减和频散累积将会影响成像的振幅和相位,进而降低成像分辨率,致使地下构造定位不准.本文基于标准线性体模型的黏弹性波动方程,提出一种衰减补偿的黏弹性最小二乘逆时偏移成像方法,推导了准确的黏弹性最小二乘理论框架下的偏移算子、反偏移算子、步长,并给出了算法的数值实现流程.通过对层状模型和Marmousi模型进行测试,得到本文提出的黏弹性波最小二乘逆时偏移方法的成像结果和收敛速度,并与其他方法的测试效果进行了对比.结果表明,黏弹性最小二乘逆时偏移方法能够准确地归位包含强衰减的多分量地震记录,能够有效避免弹性波最小二乘逆时偏移成像错误的反射层位置,压制成像噪声,收敛速度快,成像精度高.     

海冰非均匀条件下大气环流模式和海洋环流模式耦合的实现

通过采用一种较好地考虑了海气通量交换对下垫面特征非线性依赖的通量计算方案,利用较成熟的逐日通量距平耦合方案实现了次网格尺度海冰非均匀条件下一个全球大气环流模式和一个全球海洋环流模式在高纬地区的耦合,并完成了一个50年的长期积分。模拟结果表明,通量计算方案中水道的作用是合理的,耦合模式能较好地模拟出北半球高纬海冰的主要地理分布特征,但夏季北冰洋内部靠近欧亚大陆部分边缘海区海冰密集度偏大。     

塔北、塔中地区碳酸盐岩及其原油中的二苯并噻吩类化合物

塔里木盆地塔北、塔中地区四口井(库南1井,轮南46井,塔中12井和塔参1井)寒武-奥陶系12个碳酸盐岩烃源岩(泥灰岩,泥质灰岩和云岩)岩芯样品抽提物芳烃组分中的二苯并噻吩类化合物组成和丰度变化特征可分为三种类型:Ⅰ二苯并噻吩、甲基二苯并噻吩型;Ⅱ二苯并噻吩、甲基二苯并噻吩、二甲基二苯并噻吩+三甲基芴混合物型;Ⅲ二苯并噻吩、甲基二苯并噻吩、二甲基二苯并噻吩和三甲基二苯并噻吩型。研究的塔北、塔中隆起11个海相油二苯并噻吩类化合物分布类型均为二苯并噻吩、甲基二苯并噻吩、二甲基二苯并噻吩和三甲基二苯并噻吩型,和海相烃源岩二苯并噻吩类化合物第Ⅲ种分布类型完全相同。据此推断:塔北隆起8个海相油可能主要来源于塔北轮南地区下奥陶统;塔中隆起3个海相油可能来源于塔中地区中-上奥陶统。     

褶皱构造的曲率分析及其裂缝估算--以江汉盆地王场褶皱为例

曲面曲率的几何意义可以反映曲面一点邻近区域的具体形态和变形特征,因而可以用来研究褶皱层面的变形特征、构造样式及其可展性.褶皱面裂缝发育主要为一组张裂系,张裂隙走向与最大主曲率的方向垂直.褶皱构造的主曲率在一定程度上可以反映与褶皱有关的张裂隙的发育情况.对于圆柱状褶皱可以估算与褶皱有关的裂缝孔隙度和渗透率.以王场褶皱构造面为例,通过计算二次趋势面的特征值和特征向量得到褶皱构造面主曲率的大小和方向,确定该褶皱为圆柱状褶皱,并对与之相关的裂缝孔隙度和渗透率进行了估算.     

二维射线跟踪技术中的掩星面坐标系

在大气球对称假设下,可以从弯曲角观测资料通过Abel积分变换反演得到折射率廓线。在低层大气中,折射率水平不均匀性可能导致Abel反演在近地面的误差达到10%。避免该误差的一种有效方法是不通过Abel反演而采用非对称模式下二维射线跟踪算子直接将弯曲角同化到数值天气预报模式中。文章介绍了非对称模式下二维射线跟踪技术的背景、原理和数学模型,讨论了掩星面坐标系与地心惯性坐标系之间的转换关系。     

铧厂沟金矿床区域韧性剪切带特征

首次对铧厂沟金矿床区域韧性剪切带进行了较为系统的研究。根据野外地质调查和室内显微构造分析 ,区内发育一条较大区域韧性剪切带 ,无论沿走向还是顺倾向均呈舒缓波状 ,强变形带和弱变形域呈镶嵌形式。区域韧性剪切带经历了右行—左行—右行多期 (次 )活动 ,剪切方位也多次变化 ;早期形成温度约 5 0 0℃ ,以右行剪切为主 ,古应力值大于 0 .0 75GPa。控矿韧脆性剪切带是区域韧性剪切带演化的产物 ,最终形成脆性断裂。区域韧性剪切作用控制矿床、矿带的分布 ,并使部分金从矿源层分溢出来 ,产生第一阶段金的富集。次级韧脆性剪切带 (控矿剪切带 )控制富矿体的分布。因此 ,铧厂沟金矿床可称为韧性剪切带型金矿床。另外 ,中 -下泥盆统三河口群第一岩段第一岩层 (D1 -2 SH1a)部分原岩有明显海底热水同沉积特征 ;在矿床之西万家山—张家山应注意寻找硅化石英粗糜棱岩型金矿石。     

气候变化对全新世以来中国东部政治、经济和社会发展影响的初步研究

将社会科学(文史考古研究)与自然科学(气候变化)结合,分析了黄土高原社会经济发展较长江三角洲地区为早的原因,提出这可能是由于全新世大暖期的气候变化对不同地理区域产生不同的影响,最后并分析了中国历史上的3次移民潮的原因及其对长江三角洲社会经济发展的影响。