双侧向测井有限元正演算法的GPU并行实现

来源 :2014年中国地球科学联合学术年会 | 被引量 : 0次 | 上传用户:beargtg
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文根据非常规油气中的页岩气储层的地质与物性特点,建立了页岩储层的地层模型;对双侧向有限元正演算法进行了深入研究,发现其中求解线性方程组的共辘梯度算法比较耗时,因为它涉及大型矩阵与矩阵、矩阵与向量相乘,并且也发现大型稀疏矩阵的存储是很耗费资源的。利用压缩稀疏行CCSR存储技术对大型稀疏矩阵进行存储;根据矩阵相乘的特点,对其进行GPU并行优化,实现了求解线性方程组的并行共辘梯度算法;与串行算法进行对比,发现随着迭代次数的增多,串行算法所需的时间急剧增加,而并行算法所需时间增长的并不快,迭代十次以上,并行算法至少比串行算法节约50%的时间。
其他文献
大庆长垣油田高密度三维地震采集是目前国内面元最小、道密度最高采集之一,获得了高品质原始地震资料,经过与原来采集的老资料对比,基于DSU3数字检波器的高密度三维数据能够较清晰的记录信号的波场特征,原始资料的频宽比老资料频率高,且假频成份少。最终高密度成果剖面显示针对开发尺度下的小断层和砂体形态细节刻画方面都优于老资料,而且与井的符合率大大的提高,从而体现了高密度三维地震勘探技术具有巨大优势和广阔的应
粘弹性地层的吸收衰减效应使地震波在传播过程中振幅发生衰减、相位发生畸变,严重降低了地震资料的分辨率.随着地震勘探精度要求的逐步提高,地震资料高分辨率处理变得越来越迫切.综合考虑激发因素和子波的选取原则及检波器的频幅响应,本文所给的地震子波频谱表达式能较好的吻合实际情况;利用本文所提出的改进的频谱多项式拟合法或改进的谱比法,先求取地层对地震波的平均吸收系数,并进行多项式拟合后,可反演稳定的层Q值;利
实际介质的非完全弹性特征促使基于介质理论的地震波场传播特征的研究.裂缝及孔隙发育区域通常被流体充填,从而地球介质表现出强各向异性和粘弹性特征.基于各向异性介质和粘弹性介质的地震波传播特征的研究也为油气储层勘探开发提供了理论基础.粘弹性介质波场传播系统在逆时反向传播方面表现出不稳定状态,常规适用于声波波动方程的最后时刻波场反传的方法不再适应于粘弹性介质系统。同时相同时间步长下粘弹性介质系统需要更多的
本文对多分量信号应用VMD原理进行分解和对地震资料应用VMD方法进行压噪处理,并与EMD方法进行了对比,测试结果表明: VMD能够更准确的分解信号,只要参数选择合适,不会出现模态混叠现象;噪声测试表明VMD方法比EMD抗噪性更好,且去噪效果比较令人满意;由于VMD算法是在频率域进行的,其运算效率比EMD高很多。VMD是一种非常类似于EMD的强大的信号分析工具,它有着牢固的数学理论基础,这是EMD所
Russell流体因子是判识储层含流体特征的敏感指示参数。应用贝叶斯框架进行基于随机地震反演的Russell流体因子直接估算,是开展流体识别的可靠途径。本文利用序贯高斯模拟方法得到基于地质统计学的先验信息,结合确定性反演方法构建似然函数,然后引入Metropolis算法对后验概率密度进行抽样,得到反演问题的解。相比确定性反演方法,随机反演结果的分辨率比较高。并且由于序贯高斯模拟采用新的实现方式,计
塔里木盆地为中国的最大含油气内陆盆地,油气的资源潜力巨大.塔中地区则是塔里木盆地油气勘探的重点区块.而卡1区块又是中石化在塔里木盆地新登记区块中资源丰度最高的区块之一.进行卡1区块鹰山组地震沉积相研究,将为勘探区带的优选和评价提供基础资料和科学依据.本文旨在运用层序地层学、沉积学、石油地质学等理论和方法技术,以塔中隆起西北倾末端卡1三维区块为研究区,以奥陶系鹰山组为研究对象;在此背景下,充分利用露
在野外数据采集过程中,由于地形条件或者经济成本的限制,地震勘探数据沿空间方向上通常进行不规则采样,从而导致采集到的地震数据不规则、不完整,出现空间假频,使得这种稀疏分布的地震数据难以满足后续处理的要求,因此需要发展较好的插值方法,重建出缺失的地震道.另外一方面,如果有适当的插值方法,可以降低采集成本,提高施工效率.为此,本文提出直接对频率切片进行处理,也即先采用傅立叶变换,将重建数据t-x-y域变
地震数据插值和重建在地震数据优化和处理中起到了很重要的作用,不进行处理,会使覆盖次数较低的不规则或缺失数据存在成像噪声。本文讨论了一种有效的基于傅里叶变换的五维重建方法,通过在地震数据全空间内插值,可以比较接近实际地预测振幅信息,即使大间距采集的数据也可以得到好的重建结果。
油气工业界采集的地震资料大部分为三维资料,频率域全波形反演的提出使得全波形反演在3D实际资料处理上具有了可行性.全波形反演对计算资源的要求是制约其发展的瓶颈之一.近年来随着计算机技术的发展,3D全波形反演实际资料的应用开始陆续出现.全波形反演包含两个过程即,正演过程与反演过程.三维情况下仅仅使用时间域算法进行全波形的正反演,其反演过程需要巨大的存储空间,耗费较长的CPU时间;而仅仅使用频率域算法,
弹性波的纵横波分离一直是矢量波成像的一大难题,对于各向同性介质,一般在空间域使用亥姆霍兹分解,可得到较好的结果,但在各向异性介质下,由于介质的各向异性差异,导致波震动与传播方向发生了偏离,纵波不再沿着波的传播方向,横波不再垂直波的传播方向,所以如果使用亥姆霍兹分解对其进行分离将会有大量残留导致最终成像结果中出现干扰,这就会对精细成像和小构造成像带来困扰.本文将使用在频率波数域求解Christoff