基于GPU并行加速的叠前逆时偏移方法

为了提高复杂地下介质的成像精度和偏移算法的计算效率,提出可高效对地下复杂构造进行准确成像的GPU加速叠前逆时偏移方法.该方法采用双程声波方程进行波场延拓,突破倾角限制,借助于高阶有限差分方法实现叠前逆时偏移成像;利用GPU(Graphic Processing Unit)并行加速技术对波场延拓和成像进行计算,相比于传统算法,其计算效率有较大提高,可以解决叠前逆时偏移算法计算量过大问题;在获取波场信息过程中,也采用随机边界条件,实施以计算换存储策略,解决逆时偏移计算中的海量存储问题.模型测试结果表明,该方法能够高效和高精度地对地下复杂地质体成像.     

拉梅参数直接反演方法在东海N构造致密储层烃类检测中的应用

东海N构造主要目的层为强水动力环境下发育的三角洲平原分流河道砂体。平面砂层分布不连续,横向非均质性强。受埋深压实和成岩作用等影响,储层呈低孔渗特征,岩石物理规律叠置严重。此外,深部地震数据存在缺乏大角度入射信息等问题,落实研究区致密储层流体分布范围对于勘探开发设计部署具有重要意义。本文引入一种基于拉梅参数直接反演的地震流体描述方法,通过实测井数据的岩石物理定性和定量判析,优选流体敏感参数,进而结合拉梅参数的两项AVO模型参数化方程,从叠前道集中分别提取拉梅参数的AVO属性体。然后利用有色反演技术直接提取流体敏感弹性信息,以指导地震流体描述。实例应用表明,该方法的烃检结果与测井解释成果匹配度高,能够有效刻画研究区致密储层流体展布规律,可为新领域油气资源发现提供重要技术支撑。      

随机介质的声散射(Ⅰ)——广义Born近似反演

运用泛函微分方法给出弱随机介质中散射问题的散射幅度Ｂｏｒｎ近似的推广形式。讨论了关于阻尼最小二乘法反演介质参数的各价相关函数的方法,并讨论具有有限支集势函数扰动的逆问题的唯一性,其结果对运用Ｎｅｗｔｏｎ法等数值方法有实际意义。     

煤变质作用对煤层气赋存和富集的控制--以沁水盆地为例

通过对国内外有关煤变质作用与煤层气形成、赋存和富集关系研究成果的系统总结，阐述了煤变质作用类型及煤变质作用程度与煤层气的生成、吸附量、煤储层孔隙—裂隙系统形成发育及煤层特性的相互关系，认为煤变质作用对煤层气赋存富集的影响集中体现在它对煤储层孔隙一裂隙系统形成发育过程的控制上。并以沁水盆地为例，探讨了高阶煤煤变质作用对煤层气赋存富集的控制作用，提出了对沁水盆地高阶煤煤层气勘探开发研究的几点思考和建议。     

计算视复电阻率的新的近似公式

本文推导并评价了一个新的计算视复电阻率的近似公式,其视振幅表达式为|ρs（ω）|≈ρs（ω）|ρs（ω）=ρs（ω）|,视相位表达式为φs（ω）≈sum from i=1 to n（Bi（ω）φi（ω））,这里Bi（ω）是按各介质复电阻率振幅计算的稀释系数。此近似公式比迄今为止已发表过的同类近似公式都精确,可在频率域激电法数值模拟计算视复电阻率时应用。     

旋转扰动流体中壁面湍流不稳定及其拟序结构研究Ⅱ：数值模式

对第一部分的实验结果进行了讨论,指出局地最大、最小涡度的存在与所谓的正压不稳定和大气中经常观测到的滚动涡有关。速度廓线中的拐点和涡度极大值点对应不稳定发展的位置。另外,为了改进湍流的模拟,根据“准正则”近似,发展了一个三阶闭合模式。这个模型已被证明可以描述与湍流过程有关的流体动力学变量的趋势。     

变参考慢度Born近似傅氏偏移

针对常规Born近似傅氏偏移方法对于剧烈横向变速介质不能精确成像的状况而提出了变参考慢度Born近似傅氏偏移，理论上解决了任意速度变化地质模型的偏移成像问题。此外，为进一步提高复杂地层的成像精度和波场延拓算子的稳定性，对散射波场的计算公式作了改进。将改进的方法运用于盐丘模型的正演和偏移试验，并与常规Born近似偏移方法相比较，可明显看出变参考慢度Born近似傅氏偏移方法在效果上要优于后者，其处理速度横向变化的能力大大增强。     

时变轴力下结构P-Δ效应高精度分析方法

传统的P-Δ效应分析方法多为基于离散结构模型的方法,对于时变轴力情形,通常难以保证计算精度或计算效率方面的需求。为更精确地反映P-Δ效应对结构动态响应的影响,本文应用微分求积(differential quadrature, DQ)原理,提出一种高精度、无需迭代计算的结构动力P-Δ效应分析方法。以质量和刚度物理特性连续分布的大高宽比结构体系为研究背景,建立横向和轴向均任意分布的动态荷载共同作用下的结构运动偏微分方程。利用DQ原理,基于控制微分方程的强形式进行空间离散,并利用DQ权系数矩阵修正法处理边界条件。结合无条件稳定的Newmark常平均加速度法,建立包含P-Δ效应影响的结构动态响应稳定分析方法。将本文方法计算得到的动力响应与通用有限元软件的分析结果进行对比,表明在计算连续体结构P-Δ效应时,本文方法具有很高的计算精度,无需划分单元,只需使用9～10个Chebyshev-Gauss-Lobatto节点(CGL节点)即可得到非常精确的动态响应分析结果,而传统P-Δ效应分析方法需要划分近100个有限单元才能达到同样的精度。本文方法还能清晰地揭示时不变轴力和时变轴力改变结构横向振动特性的...     

一个双波地形重力波拖曳参数化方案

当地形次尺度强迫的作用与显式的经典动力作用效应相当时,地形重力波拖曳力对于环流的维持,以及动量和热量通量输送的动力效应变得十分显著。这种地形次尺度拖曳作用项可通过参数化的方法,在动力方程中加入额外的小项而引入数值模式。目前成熟的地形重力波拖曳参数化方法,如第1代基于线性单波理论的参数化方案;以及侧重考虑了临界层作用等因素对拖曳力的额外贡献的第2代参数化方案,都无法有效表达风速垂直变化引起的波动应力随高度变化的特征。基于上述考虑,本文给出了一个双波参数化方案用于计算地形重力波拖曳中由线性自由传播重力波造成的波动应力的垂直分布。通过二阶WKB近似,它对由风速垂直变化引起的对波动应力的选择性临界层吸收过程和经典的临界层吸收过程做了显式处理;而在不发生临界层吸收现象的地区,则用两个单波同时在垂直方向上进行应力的传播,并利用波饱和标准进行应力耗散。进一步地在真实地形(以大别山地区为个例)条件下的测试结果表明,通过在不同理想风速廓线以及北半球冬季中纬度纬向平均风廓线下对波动应力垂直分布的计算,证明该方案确实能有效地给出应力随高度变化的特征。