基于Remez迭代算法求解差分系数的双相介质自适应有限差分正演模拟

利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性.     

用柱状成层各向异性介质的并矢Green函数模拟偏心对多分量感应测井响应的影响

采用偏心状态下柱状成层各向异性(横向各向同性)介质中并矢Green函数的解析表达式高效模拟多分量感应测井仪器在井眼中偏心时的响应.为提高精度,在模拟时考虑到了金属心轴、绝缘保护层的存在以及各分量线圈系的具体形状.数值模拟结果表明,当井眼内钻井液电导率相对较高、地层电导率相对较低时,偏心对仪器响应的影响较大,尤其是对短线圈距线圈系的影响更为明显,必须进行偏心效应校正.当钻井液电导率相对较低时,偏心对线圈系响应的影响可忽略不计.对位于相对低电导率井眼中的线圈系而言仪器方位角的影响可忽略不计,而当线圈系位于相对高电导率井眼中时仪器方位角的影响极为明显.当仪器偏心率较小时线圈系的响应随仪器方位角的变化较小,仪器偏心率越大线圈系的响应随仪器方位角的变化越明显.     

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.     

基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移

近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制(η0)和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明:①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略.     

基于全局弱式无单元法直流电阻率正演模拟

全局弱式无单元法是在有限单元法基础上发展起来的一种数值模拟方法,它采用局部支持域内的节点信息来构造形函数实现局部精确逼近,摆脱了单元,仅依赖于节点信息,具有预处理简单、模拟精度高、灵活性强的特点,适用于复杂地电条件下直流电阻率正演模拟.本文采用RPIM构造直流电阻率全局弱式无单元法形函数,利用RPIM形函数推导了直流电阻率全局弱式无单元法方程.然后,编制了直流电阻率全局弱式无单元法正演模拟Fortran程序,利用该程序对典型的地电模型进行了正演模拟,并将正演结果与基于线性插值的FEM正演结果及解析解进行对比,结果表明采用RPIM形函数的全局弱式无单元法用于直流电阻率正演模拟的正确性及有效性,且在同等条件下,全局弱式无单元法模拟精度高于矩形剖分的FEM,更有利于指导电法勘探的数据解译;利用该程序对复杂地电模型进行了正演模拟,结果表明全局弱式无单元法对复杂地电模型模拟效果良好,适应性强,灵活性高,可任意加密节点提高模拟精度.     

大地电磁三维正演并行算法研究

大地电磁法三维正演算法计算量大,用传统的串行程序计算相当耗时。而三维正演是逐个频率按顺序计算的,并行性好,适合并行运算。结合MPI自身的优越性,在深入分析大地电磁三维正演串行程序实现流程的基础上,确定了并行计算的思路,实现了三维正演的并行计算。通过三个理论模型对实现的三维正演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维正演并行程序运行结果正确,效率提高明显。此思路可为解决其它地球物理超大计算量问题所借鉴。     

正二十面体六边形全球离散格网编码运算

全球离散格网系统是指把地球表面按照一定规则离散分割成多分辨率层次结构的格网单元,广泛应用于海量多源空间数据的组织、管理和分析中。六边形全球离散格网具有优良的几何特性,非常适合于空间数据的处理,如何进一步提高六边形全球离散格网编码运算的效率仍是当前研究的重点。本文采用正二十面体施耐德投影四孔径六边形全球离散格网模型,基于六边形三轴坐标与编码的二进制数的对应关系构建四孔六边形的基础编码结构,将二十面体划分为32个基础六边形,并将之分为3种基础六边形剖分瓦片,在每个六边形剖分瓦片采用基础编码结构进行编码,建立了四孔六边形全球离散格网编码,同时设计了并实现了四孔六边形编码与六边形三轴坐标之间的快速转换,基于此构建了一种高效的四孔六边形全球离散格网编码运算方案,包括编码的算数运算、空间拓扑运算和邻域检索运算及跨面运算。与现有的六边形全球离散格网编码运算方案相比,本文的方案进一步提高了编码算数运算、空间拓扑运算和邻域检索运算的效率,编码加法运算是HLQT的2～3倍,邻域检索运算分别是HLQT的3～5倍和H3的2～3倍,且受格网编码层次的影响较小,编码的跨面邻域检索运算时间略高于面内的运算,可以为全...     

基于二次场的海洋可控源电磁法三维有限元正演

根据库伦规范势的定义,推导出关于磁矢量势和电标量势的偏微分方程,为了克服由电流源引起的奇异性和数值模拟计算困难,将电磁总场分解为一次场和二次场,一次场由基于Schelkunoff势函数的一维正演算法得到,二次场由有限元法计算得到,实现了海洋可控源电磁法三维有限元正演算法。通过一维数值模拟实例,验证该算法的计算精度。然后,利用该算法对带海底地形的三层储层模型进行正演,分析了海底地形对海洋控源电磁场各分量产生的影响。      

基于Ruger近似式预测裂缝型储层的可行性研究

叠前纵波方位各向异性检测裂缝是目前应用最为广泛的方法,该方法主要是利用AVAZ地震资料椭圆拟合的长轴与短轴信息来评价裂缝,该方法在具体应用时,认为裂缝型储层AVAZ响应特征的变化仅由各向异性参数决定,而没有考虑非各向异性参数的影响,从而引起裂缝评价的多解性。为提高该技术的裂缝评价精度,从模型分析入手,将模型中各向异性参数及岩性参数设置为概率密度分布函数,然后采用Monte Carlo随机方法进行叠前AVAZ正演模拟。首先,对各向异性参数对椭圆拟合的影响进行简单分析,得出各向异性参数γ对椭圆扁率B/A与各向异性因子B影响最大、δ次之、ε最小。然后,进一步重点分析不同标准差的速度与密度模型的B/A与B响应特征,得出地层的速度与密度的变化在较小范围内时,该技术评价裂缝才具有一定的可行性,其中,地层纵波和横波速度的变化（尤其是纵波速度）对叠前AVAZ响应影响较大,地层密度几乎没有影响,故当横向岩性变化较大时,该方法在评价裂缝型储层时的结果具有一定的不确定性。      

中条山地区涑水杂岩新太古代烟庄正长花岗岩年龄及成因:对华北克拉通地壳演化的制约

太古宙末期钾质花岗岩的广泛发育是陆壳成熟和稳定化的重要标志,对了解早期陆壳的形成与演化具有重要的意义.发育于华北克拉通南缘中条山地区涑水杂岩中的烟庄正长花岗岩的形成年龄和成因还没有被很好地限定,构造背景还存在争议.对烟庄花岗岩进行了锆石U-Pb年代学和Hf同位素以及全岩地球化学和Nd同位素研究.烟庄花岗岩锆石SHRIMP U-Pb年龄为2 515±7 Ma.岩石具有高硅(SiO₂=73.05%~74.85%)、高钾(K₂O=4.46%~5.86%)、富碱(ALK=8.32%~9.36%)、贫钙(CaO=0.55%~0.98%)、低TFeO*(0.73%~1.28%)和MgO(0.31%~0.52%)的特征,A/CNK=1.01~1.04,为弱过铝质的钾玄系列.稀土总量变化较大(ΣREE=63.80×10-6~250.02×10-6)),轻重稀土元素分异明显((La/Yb)_N=25.44~92.87),Eu异常变化较大(Eu/Eu*=0.47~0.79).岩石低Sr、Ba,富集Rb、Th、U等元素,亏损Nb、Zr、Y、Yb、Cr、Co、Ni等元素,具有较高的Rb/Sr、Rb/Ba和Sr/Yb比值以及较低Sm/Nd和Nd/Th比值,具有高分异I型花岗岩的特征.烟庄花岗岩具有0附近的全岩ε_Nd(t)值,岩浆锆石具有正的ε_Hf(t)值(2.85~3.66),两阶段Hf模式年龄为2 258~2 883 Ma,多数在2 600~2 883 Ma之间.结合其他方面研究,烟庄花岗岩具有同碰撞和后碰撞花岗岩的特征,推测为新生地壳在由挤压向伸展转换的构造背景下部分熔融所形成,可能有少量地幔物质添加.这期钾质花岗岩的形成,标志着华北克拉通太古宙末期强烈岩浆活动的结束以及稳定陆壳的形成.