弹性波方程正演的粗粒度并行算法

波动方程正演在地震资料采集、处理、解释与反演中均能发挥重要作用,但现有的基于求解地震波动方程的正演算法由于受庞大计算量的制约而难于大规模应用于工业生产.本文从二维弹性波动方程出发,研究了利用有限差分法并行求解该方程的基本思路与方法,给出了适于并行求解的计算空间划分方法与通信方案,分析了不同参数条件下并行程序的运行时间、加速比与效率.引入消息传递接口（MPI）实现了弹性波方程的并行求解,极大地提高了数值求解弹性波方程的计算效率.     

基于单元细分H-自适应有限元全张量重力梯度正演

在重力梯度各阶张量正演中,难于求取复杂地质体的解析式,转而利用已知解析解的简单规则几何形体近似剖分地质体,其会引起较大的几何拟合误差.本文首先从全张量重力梯度解析公式出发结合等参变换基于有限元分析给出了满足地下复杂地质体正演的型函数,引入体积分的误差指示器用以估计全张量重力梯度误差,提出了适用于全张量重力梯度正演的自适应计算的策略及迭代算法,数值算例证明了本文方法的正确性和有效性.     

正逆累积法和分布函数法在研制储层电性标准中的应用

针对储层电性标准研制时采用的交汇图法--单一的 "L"型划分方法存在的不直观性和人为因素的影响,将正逆累积法和分布函数法首次应用到确定电性标准图版中,并以喇嘛甸油田气顶区为例,进行了气砂(含气砂岩)电性标准的研制.误差统计表明,两种划分方法均达到了国家储量规范的要求.这使得电性标准划分更加合理、精确与直观,为储层电性标准的研制提供了新的方法与思路.     

某区块三维深水采集设计与论证

某区块主要勘探目的层位于2.5～3.5 s,埋藏较深、断裂发育、构造破碎、各年度存在闭合差.地震资料的信噪比与分辨率不满足油气勘探需求.为此,针对该区块的地震地质条件,结合照明分析以及兼顾各项采集设计要求与工作效率的条件下,重新论证并设计了三维地震观测系统,并给出了海上三维采集设计的处理思路.     

考虑影响区真实形状的竖向排水井地基固结解

针对现有竖向排水井地基固结解析理论对影响区一律采用圆形等效假定的缺陷,研究了影响区真实形状为正六边形的按梅花形布置的竖向排水井地基的固结问题。通过建立新固结方程和引入新的边界条件,并考虑土体水平渗透系数变化,得到了相应的解析解。对于土体水平渗透系数的3种变化模式,分别给出了各种模式下的特殊解。在此基础上,分析了3种模式下3个主要的无量纲参数对地基固结性状的影响,并比较了计算结果和现有理论结果。分析结果表明,影响区和涂抹扰动区范围越大,固结越慢;土体的最大与最小水平渗透系数之比越小,固结越快;在相同条件下,考虑扰动区渗透系数线性变化的模式2固结最快,而假定扰动区渗透系数不变的模式1的解与现有理论解相当接近,验证了现有竖向排水井地基固结解析理论中对影响区采用圆形等效假定的合理性     

二阶声波方程频域PML边界条件及频域变网格步长并行计算

研究了二阶声波方程频域PML边界条件和频域变网格并行计算技术。PML边界是一种较为理想的吸收边界方法,多用在求解时域应力速度方程中,但对于频域声波正演,二阶位移方程更常用。从一阶声波方程PML吸收边界条件导出频域二阶位移方程PML边界条件,模拟算例得到的频率切片、时间切片和地震记录对比都说明该边界条件吸收效果很好。频域单炮正演不同频率间是独立的,据此低频部分采用大网格计算,高频采用小网格,实现变网格步长计算技术,这是较时间域正演的一个优势,在保证模拟质量的同时,减少计算量和内存消耗。     

掘进巷道直流电法四极装置超前探测的模拟及应用研究

为研究矿井直流电法四极装置在掘进巷道超前探测中的效果,本文以矿井直流电法超前探测原理为基础,通过分析COMSOL Multiphysics点电源正演模型模拟值与理论值的误差,验证了COMSOL Multiphysics研究直流电法超前探测的可靠性与可行性;建立巷道超前探测三维地电模型,通过正演模拟,分析了巷道空腔对数据的影响,对比了二极、三极装置与四极装置的探测效果及对不同距离低阻体的响应特征;基于相似性原理构建物理模拟实验,绘制了超前探测视电阻率等值线图,研究了四极装置对迎头前方不同形状、距离异常体的实际探测效果.研究结果表明：直流电法四极装置能够对迎头前方低阻体产生异常响应,探测数据剔除巷道空腔影响后,能够精准识别出异常体的横向位置;并且相较于板状异常体,四极装置对球状异常体的超前探测效果更好．通过将直流电法四极装置超前探测应用于工程实例中,经钻探验证与探测效果一致,说明该技术具有很强的实用性和推广性．     

可控震源粘弹性波动方程有限差分模拟

在一般的地震波计算中,地球介质可以当作各向同性的完全弹性体来对待。而对地表土壤和精细观测,粘弹性介质模型比完全弹性模型更符合实际。可控震源是近年来发展起来的精确可控、安全有效的新型震源,但国内对可控震源传播动力学的研究较少。对可控震源激发的地震波在粘弹性介质中的传播,用有限差分方法进行了数值模拟,给出了Kelvin一阶粘弹性波动方程组、吸收边界条件和震源处理。提出了检波器接收信号加速度和传播距离的三次多项式指数关系式,比线性指数关系更准确。使用该公式可以预测采集系统有效接收距离,指导测线布置。      

斜井三维VSP多波射线正演方法

基于费马原理和 Snell 定理,对三维 VSP两点射线追踪问题作了研究,给出了三维逐段迭代射线追踪算法的计算公式;同时考虑射线入射到界面时能量分配问题,探讨了三维介质条件下射线追踪计算方法。对斜井VSP三维任意起伏层状介质中多波进行正演模拟,结果表明:该方法可行性和有效性较高,计算速度快,且计算精度高。      

内蒙古敖汉旗大黄花正长花岗岩锆石U-Pb年代学及地球化学特征

通过研究大黄花正长花岗岩的锆石U-Pb年代学,结合地球化学特征,探讨其形成时代、岩石成因及其构造背景.大黄花正长花岗岩的同位素测试结果为162.6±1.9 Ma,表明该岩体形成于中侏罗世晚期.该花岗岩主量元素具有高Si、富碱、低P和Ca的特征,微量元素具明显的Eu、Ba、Sr、P、Ti负异常.岩石经历了高分异演化（DI=95.1~95.88）,为高分异I型花岗岩.极低的Sr/Yb比值暗示其形成于一个非常低压的熔融环境.结合岩石地球化学、区域地质特征,认为大黄花正长花岗岩是蒙古-鄂霍次克缝合带演化的产物,其形成于碰撞后的伸展环境.