改进的散射积分算法在多参数波形反演中的应用

来源 :2015年中国地球科学联合学术年会 | 被引量 : 0次 | 上传用户:Robert_1967
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
多参数全波形反演是当前研究的热点,地震数据对不同参数的敏感性不同使得不同参数在反演过程中相互耦合,增加了反演的非线性程度。反演策略在一定程度上可以降低多参数全波形反演的非线性,然而,反演策略并不能从根本上解决不同参数间的耦合问题。Hessian矩阵的有效利用是解决多参数全波形反演非线性问题的有效途径。Liu等提出的BKFWI方法可以方便地扩展到多参数全波形反演中,通过Hessian矩阵的迭代逼近实现多参数之间的解耦。因此,本文分别基于VTI介质拟声波方程、各向同性弹性波方程、变密度声波方程,导出了不同参数所对应的核函数,并阐述如何利用BKFWI算法实现各向异性参数、纵横波速度、速度与密度的多参数联合反演。全波形反演(Full Waveform Inversion,FWI)在最优化理论框架下,通过拟合模拟数据与观测数据,使残差达到最小,从而获取具有高分辨率的地下物性参数.梯度的计算是全波形反演的核心,根据散射积分原理,在最小二乘目标泛函下,目标函数的梯度总可以表达为敏感核函数与数据残差相乘的形式.本文所采用的改进的散射积分算法无需存储核函数与Hessian矩阵即可实现各种方向与步长的计算,尤其在高斯-牛顿方向的构建中无需额外的正演计算,与目前广泛采用的伴随状态法相比,显著提高了高阶全波形反演的计算效率。另一方面,多参数全波形反演中,不同参数之间的耦合效应包含在Hessian矩阵中,截断高斯-牛顿算法可以实现对Hessian矩阵的有效利用,实现多参数的解耦,上述数值实验证明了本文的观点。
其他文献
晚古生代古亚洲洋的闭合时限及其相关块体的碰撞拼合过程一直存在很大争议,内蒙古中东部地区位于中亚造山带东南缘,对于研究中国东北诸块体碰撞拼合过程及其与华北板块拼合形成华北-蒙古统一陆块的演化历史具有关键意义[1][2]。本文对内蒙古中东部晚古生代地层进行了系统的古地磁学采样,对820块样品进行了详细的岩石磁学分析,采用交变退磁与热退磁,或二者混合退磁的方法,从380块样品中分离出了特征剩磁分量(Ch
以大别山双河地区超高压硬玉石英岩为主要研究对象,应用傅里叶变换红外光谱仪(FTIR)、电子背散射衍射(EBSD)及电子探针(EMPA)对硬玉和石英的含水量、变形组构及地震波性质进行了详细研究,同时对比研究了阿尔卑斯山超高压变质带意大利Sesia 地区硬玉-钾长石-多硅白云母-正片麻岩中的硬玉和石英的含水量和变形组构.红外光谱研究结果表明,硬玉的吸收峰主要分布在3350-3370cm-1、3452-
江南造山带是一个中-新元古代的造山带,理解江南造山带的形成过程对理解华南前寒武纪的地壳演化具有重要的意义。然而,江南造山带东段新元古代地层序列发育具有巨大差异性,高精度同位素年代学数据依然缺乏;此外,由于缺乏典型剖面精细的沉积相研究,沉积物物源仍不清楚,制约了对沉积盆地性质的认识。本文以浙北地区河上镇群和皖赣相邻地区的双桥山群为研究对象,运用地层学、沉积学、岩石学、碎屑锆石U-Pb 年代学和Hf
蛇绿岩是古大洋岩石圈的残片,并通过板块俯冲作用拼贴到大陆边缘或残存于俯冲-碰撞造山带之中,它对研究板块构造的启动时间、机制及大地构造演化及超大陆的重建具有重要意义。作为扬子克拉通目前唯一出露的太古宙TTG 陆核,黄陵穹窿核部是研究扬子克拉通早期地壳演化的重要窗口。最新的研究表明,该区古元古代存在碰撞造山变质与造山后裂解作用事件有关的重要记录。本文报道了出露于扬子克拉通黄陵穹窿北部太古宙-古元古代杂
Hf(t)(–16.7~–32.7)值,暗示它们起源于下地壳物质的部分熔融.此外,曲里早白垩世晚期低Mg#埃达克质岩石中新元古代(754–542 Ma)岩浆锆石和三叠纪(220±4 Ma,n = 9)变质锆石的存在,结合相对高的206Pb/204Pb 初始比值暗示,它们起源于经历了超高压变质作用改造的扬子克拉通基底物质的部分熔融.古元古代和太古代继承锆石在早白垩世早期低Mg#埃达克质岩石中的存在,
The Qilian–Qinling region of central China located in the northern margin of the Tibetan Plateau,where the Triassic sediments deposited may provide keys to understanding the tectonic evolution of the
随着油气勘探开发程度的逐步深入,成藏条件复杂和识别难度大的隐蔽性储层成为了石油勘探的重点.由于地震资料的带限性质,传统反演的分辨率很难达到要求,随机反演的应用变得更加广泛.地震随机反演可以将随机模拟技术与优化理论相结合,有效地利用测井数据中所包含的地质信息.将测井数据划分的岩性信息作为随机反演的约束,可以提高反演的准确性.利用傅里叶滑动(FastFourier Transform-Moving A
随着我国油气藏勘探的不断深入以及勘探工区难度日益加大,单靠地震技术已无法满足当前油气藏勘探需求,本文首先采用体积分方程法,正演模拟地震波场及电磁场数据。模拟计算过程中,采用稳定型双共轭梯度法(BCGS)求解网格剖分后推导的线性方程组;将并矢格林函数分解为关于场源上下层状介质界面深度的函数,如此体积分核的计算可引入卷积算子和相关算子,通过引入快速傅里叶变换(FFT),积分核的计算得到加速,计算复杂度
本文针对PSInSAR算法中整体基线网络的布设方法及具体信息存储算法进行研究,利用附有距离限制原则的Delaulanay网络布设方法和邻接矩阵模型对基线边信息进行存储.为验证算法可行性,本文利用河南地区PALSAR影像中提取出的PS点及预先安装布设的人工角反射器点开展整体基线网络布设及信息存储试验.
本文以OA介质反射系数近似公式为基础,选取了不同方位角的地震数据,研究页岩气储层各向异性梯度的提取方法,利用叠前地震数据对页岩储层各向异性梯度的提取方法进行了试算,验证了AVO属性提取方法的合理性和有效性。通过上述方法提取的页岩储层各向异性梯度的高值区域与测井数据中裂缝存在的位置有着非常好的对应关系,为页岩储层裂缝预测提供了有力的信息。但在实际操作中,方位角的选取显得尤为重要,需确保所选三个方位角