长江口斜压诊断模式三维流场数值模拟

建立了σ坐标系下长江口斜压诊断模式三维流场数学模型,采用k-kl二方程紊流闭合模型求解垂向涡粘系数,计算域内恒定、非均匀盐度场反映了计算域密度斜压效应。验证结果表明,斜压诊断模式的模拟成果较切合实际。模型成功复演了长江口水域由往复流向旋转流过渡的流场特性;潮波从四条汊道传入后在分汊口附近相遇形成多个会潮点,其位置随着径流的大小、潮汐的强弱、汊道的形态等因素而变化。     

基于地震叠后反演与属性分析的页岩气有利储层预测

中国南方页岩气勘查区块多为勘探新区,仅部署少量二维地震和地质钻井,且资料品质普遍较差,页岩气储层预测很难像成熟勘探区一样开展系统性研究工作,如何对现有数据进行充分挖掘是解决勘探新区页岩气储层预测的关键。为此,综合利用不依赖初始模型的叠后约束稀疏脉冲反演方法和多种频率类属性分析技术,开展页岩气有利储层预测。地震反演能够获得多种参数参与储层预测,属性分析能够从多种角度识别储层,通过多种信息综合应用和叠合分析以提高储层预测精确性和有效性,为新区页岩气勘探提供了一种新的思路。该项技术应用于南方某区块,圈定了页岩气的有利目标区,在该区部署的两口水平井均获得高产工业气流。     

微分算子展开与耦合势三维有限体积法快速计算随钻超深前视电阻率测井响应非线性Born逼近

本文将非均质各向异性地层中电磁场耦合势Helmhotz方程分解为各向异性背景场方程与散射场方程,并应用算子展开技术推导出耦合势Born级数解,建立一套三维非均质各向异性地层中随钻超深前视电阻率测井响应各阶散射电磁场与非线性Born逼近的有效计算方法.首先,利用传输线法与二维插值技术确定背景电磁场空间分布,得到零阶Born逼近解,然后根据电磁场耦合势的Born级数解推导出各阶散射电磁场方程与递推关系,通过三维有限体积法实现各阶散射场耦合势方程的离散,在此基础上结合MKL PARDISO并行技术,通过递推方式逐步计算各阶散射电磁场,并根据Born级数部分和确定不同阶次的非线性Born逼近.最后利用数值模拟结果对正演算法加以检验,并深入研究分析各阶散射场的变化特征与Born级数的收敛性质.

     

D-T2二维核磁共振测井反演方法对比与适用性研究(英文)

D-T₂二维核磁共振测井技术能够直观地区分不同类型、不同赋存状态的孔隙流体,在油气勘探中应用广泛。目前,,国内外学者提出了很多二维核磁共振反演方法,但大多止步于数值模拟,需要研究不同反演方法针对不同采集序列的测井数据的适用性,优选反演方法,为二维核磁共振测井数据处理的软件模块开发提供参考。本文者,对此问题,基于二维核磁共振测井原理,详细研究了TSVD、BRD、LM-模平滑和TIST-L1正则化反演方法,分别对油水模型、气水模型数值模拟得到的CPMG和DE序列的回波信号,进行反演结果精度和计算效率的对比研究。通过数值模拟,对比分析了TSVD、BRD、LM-模平滑和TIST-L1正则化四种反演方法的反演结果和计算耗时,并优选反演方法对MR Scanner仪器SP模式测井数据处理。结果表明,对于油水和气水两种模型的CPMG和DE序列回波信号,TIST-L1正则化方法精度最高,LM-模平滑方法略差,但耗时较少,更适合测井数据处理,该方法针对A25井的处理结果与Geolog软件相符,能够满足数据处理的需要。     

频率域全波形反演方法研究进展

全波形反演是一种利用地震波传播的动力学特征来获取地下介质物性参数的反演方法,可为揭示地下精细结构提供重要依据。本文以弹性波方程作为数学模型来模拟地震波传播规律并进行相应的反演方法研究。为提高计算效率与反演结果的准确性,可将近似解析离散化（NAD）算子用于频率域弹性波方程的正演模拟。本文在频率域NAD离散的基础上推导阻抗矩阵的稀疏分块结构与反演目标函数对模型参数的梯度计算公式,由此建立基于NAD算子的频率域弹性波全波形反演方法。为验证该方法的有效性,文中通过数值实验对多种典型介质模型进行反演计算,均得到了理想的反演结果。

     

地面瞬变电磁勘探中源外地形的影响

起伏地形为中国矿产资源的发现和勘查造成了极大困难,为了避免复杂的地形效应,在地面瞬变电磁勘探中尽量将发射源与接收点布置在地形平缓处,但源外起伏地形的影响仍然存在。目前,对地形影响的分析大多从测道曲线图或者视电阻率断面图进行定性判断,较少从整个观测时间段内量化评估地形影响的程度。基于有限元法实现了带地形条件下矩形回线源的瞬变电磁三维正演,在源外构造了山脊、山谷地形,将地形与发射源的距离、大地电导率、异常体电导率作为变化参数,给出了整个观测时间段内在不同接收点处地形造成的相对误差分布,从而分析了源外地形对电磁响应的影响规律。结果表明：源外地形对电磁响应的影响主要集中在某个时间段内,随着与发射源距离增大,地形作用减弱,受影响区域也推移到较晚期时间段上;随着大地电导率增大,地形主要影响区域会推移到较晚期时间上,但受影响程度基本一致;相同几何尺寸的山谷地形比山脊地形对电磁响应的影响更大;地形、异常体对电磁响应的影响可以分离,如果异常体产生的二次场比较强,地形效应的影响相对较小,可以忽略地形影响。研究结果为地面瞬变电磁勘探中地形效应的识别与后期处理解释提供了理论依据。     

地震孕育、发生、发展动态过程的三维有限元数值模拟

为探讨用数值方法模拟地震孕育、发生、发展动态过程的途径,从线性流变体介质内制约质点运动的运动方程出发,导出了模拟地震孕育、发生、发展动态过程的三维有限元方程及程序。还给出了模拟进程中地震孕育、发生、发展过程的约束条件,使得可以用同一个程序完整地模拟地震孕育、发生、发展的全过程,为地震过程本质的认识以及物理预报的实现提供了一个潜在的新手段。     

双侧向电阻率测井反演算法研究

双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难．本文利用变形玻恩（Born）迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演．在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林（Green）函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量．通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小．     

青藏高原及邻近地区地壳与上地幔剪切波三维速度结构

本文利用30个基准台所记录的238条长周期面波资料,经过适配滤波和分格频散反演,得到中国大陆及邻区147个分格10—105s的纯路径频散,进而反演出青藏高原及邻近地区深至170km的剪切波三维速度结构.研究表明,青藏高原中西部地区和东部地区的地壳平均厚度分别为70±7km和65±7km,地壳平均剪切波速度分别为3.55和3.62km/s,上地幔顶盖平均速度分别为4.63和4.61km/s; 岩石层厚度均为120±10km;东部地区下地壳内30—40km深度处普遍存在低速层;青藏高原及其东侧的上地幔低速层内有横贯东西且明显向上隆起的低速腔.滇西缅北地区的地壳厚45±5km,上地壳及下地壳内都有低速层;上地幔顶盖的速度为4.42km/s,比青藏高原本体及恒河平原都低.恒河平原地壳厚34±2km,速度平均为3.45km/s;上地幔顶盖厚86±10km,速度平均为4.63km/s,顶盖内55—83km深处有一个低速夹层.     

基于高精度DEM和RS技术的滑坡环境指标参数快速提取方法——以陕西省宝鸡市金台区长乐塬特大型滑坡为例

DEM和RS技术是研究滑坡地质灾害的重要资料和手段。近年来,随着高空间分辨率遥感卫星和高精度雷达卫星的上天,可以获取现时性高精度的DEM,使滑坡地质灾害的研究由二维向三维提升。利用IRS-P5数据生成的5m精度的DEM,借鉴GOOGLE的三维可视性原理,将其和高空间分辨率QuickBird(0.61m)数据叠置到数字地球之上,制作成三维可视性图像,进行滑坡环境指标参数提取方法研究。研究结果表明,该方法可直接读取滑坡环境指标的三维参数,具有客观、准确、快速的特点,可为滑坡灾害评估和区域地质灾害危险性评价提供定量化资料。