收藏本站
《吉林大学》 2013年
收藏 | 手机打开
二维码
手机客户端打开本文

基于Extrapolation Tikhonov正则化算法的重力数据及梯度多分量数据的3D反演方法研究

刘银萍  
【摘要】:随着重力及梯度观测数据类型的增加、硬件设施的发展、数学方法的不断进步,重力及梯度数据观测区域日趋广泛、解释技术及方法不断进步,应用领域不断拓宽,重力及梯度张量数据3D反演方法的研究也日益紧迫。反演方法的改进、计算精度的提高、计算时间及计算量的减小、可反演数据类型的丰富成为人们研究的热点。 3D反演方法分为空间域和频率域反演方法。空间域反演方法主要由两类构成,第一类是利用丰富的地质先验信息约束解,在这一分类中具有代表性的方法有种子方法,交互式反演方法等;第二类是尽量少的使用复杂地质信息的约束最小化反演方法,在这类方法中最常用的约束条件包括深度加权约束,平滑约束及上下限约束等。本文重点研究第二类反演方法。除多数学者致力研究的空间域反演方法外,以计算速度快,硬件要求较低为特点的频率域反演在特定的研究领域依然占有重要的地位,本文就这一方面也做了部分研究。 空间域反演方法的计算精度、计算时间、计算量、反演数据类型及反演约束条件等问题是近年来人们研究的焦点,本文就这些问题所做的研究工作如下。 1.深入讨论了Sorokin、Haaz、Nagy及Okabe和Steiner的棱柱体重力模型正演方法,并根据这四个正演方法计算同一棱柱体模型引起的重力异常。计算结果表明,虽然四个计算公式的最后一项不同,但是它们计算的重力异常之间的误差小于10-10,所以,四种方法的计算误差对反演结果没有明显的影响,反演中可忽略这些计算误差,选用其中任意一种正演方法进行反演计算都可满足误差水平,本文选用Nagy的正演方法计算棱柱体正演模型引起的重力值。 2.详细推导计算重力梯度九分量的计算公式,讨论了梯度数据的性质及不同分量之间的对称关系,并建立正演模型,计算此模型的梯度九分量数据,计算结果表明,梯度九分量存在很好的对称性,同时,计算结果也验证了梯度场是无迹的。 3.本文详细分析了反问题病态解的诱因,即解的存在性、唯一性和稳定性三者不能同时满足时,反演解是病态的,并给出了解决病态线性反演的方法。针对线性病态反问题,经典和简单的方法是利用Tikhonov正则化算法解决重力数据反问题中核函数的奇异性问题。为进一步减小因正则化参数的引入带到解中的误差,后人发展的Extrapolation Tikhonov正则化算法成为本文研究反问题的主要工具。 4.深入研究基于Extrapolation Tikhonov正则化算法的重力数据约束反演方法。通过比较基于Tikhonov正则化算法反演结果及基于ExtrapolationTikhonov正则化算法反演结果可知,后者的观测数据与预测数据拟合误差更小。由于正则化算法中参数选择直接关系到反演结果的好坏,文中选用平衡原则、单调误差原则及离散原则三个选择原则综合选择正则化参数n和。 5.在重力及其梯度数据反演中,核函数呈现随深度增加而快速衰减的特性,如果直接进行反演计算,则反演结果存在严重的趋肤效应。为此,前人提出一种可消除这种衰减的方法,即利用与核函数随深度变化具有相反特性的深度加权函数作用与核函数,使核函数沿垂直方向对所有密度块体的作用相等,从而消除重力数据反演中的趋肤效应。本文在深入研究前人深度加权约束函数的基础上,改进了深度加权函数。将前人加权函数、作者改进的加权函数及不加权的核函数应用到不同埋深的单个立方体模型的密度分布反演中。结果表明,不加权核函数反演的密度分布存在严重的趋肤效应,前人加权函数对于立方体的底部分辨率比改进后的加权函数要差,改进后的加权函数对于埋深不同的异常体都有较好的反演效果。 6.密度上下约束函数的引进,对超出先验信息范围的密度差起到约束作用,将反演结果转换到满足先验信息的数值范围内,文中详细分析了上下限约束函数中各参数的选择机制,并将其应用到重力及梯度多分量数据的反演计算中。 7.本文将基于Extrapolation Tikhonov正则化算法及Tikhonov正则化算法的反演方法应用于两组组合模型,就它们的计算精度及计算时间进行对比,结果表明,两个方法都可较准确的反演出异常体密度分布特征。前者能在相同的误差要求下达到更好的拟合误差水平,但是需消耗更多的计算时间。 8.详细讨论了由不同梯度分量及分量组合构成核函数时包含的信息量的计算方法及判断依据。利用奇异值分解方法分解带有深度约束等约束条件的核函数特征值,并将核函数的特征值利用最大特征值规则化,利用规则化之后的特征值谱研究不同的核函数含有的信息量。前人研究结果表明,核函数的特征值越大含有的信息量越多。将不同分量及不同分量组合成的核函数分解后发现单个分量的特征值谱小于分量组合时的特征值谱,因此梯度分量组合反演密度模型时,可获得更多垂向和水平方向的细节。 9.利用基于Extrapolation Tikhonov正则化算法的3D反演方法反演重力梯度多分量数据不同组合方式时的密度模型。结果表明单独反演不同的梯度分量时,Vxy及Vzx在水平及垂向的分辨率较理想,Vxx是各梯度分量数据反演结果中效果最差的,这与数值分析结果相一致。少数梯度分量数据组合反演结果与单个梯度分量作为核函数的反演结果相比,少数梯度分量组合的核函数有更大的特征值谱,因此,可提供更多关于异常源的信息。计算结果表明Vzz|Vzx|Vzy组合的反演结果优于其它的单分量数据及两分量数据反演结果。随着反演中梯度分量的增加,是否分量越多,反演结果越精确成为关注的焦点问题之一,为说明此问题,文中做了五组3个、4个及5个不同分量组合的反演计算,结果表明,这五组反演计算结果相近,由此可见,当参与反演的数据可同时获得较好水平及垂向分辨率时,再继续增加梯度分量数据并不能进一步改善反演结果,只能增加计算量及计算时间,带来时间和空间的耗损。 在空间域3D反演快速发展的同时,频率域重力数据3D反演的研究工作也在不断的推进,其在重磁数据解释方法中占有重要地位,就这一方面本文的研究工作包括以下几点。 1.详细论述了Oldenburg-Parker法计算重磁数据的方法原理,并利用此方法计算余弦球冠模型引起的重力值,为后文的频率域反演提供试算数据。 2.深入研究了重力数据的频率域反演方法,此频率域反演方法的主要思想如下:当等效密度层厚度与待反演的构成异常体的棱柱体水平尺寸满足一定的几何关系时,视密度与异常源深度顶部埋深之间满足拟线性关系,此时可在频率域通过线性反演计算异常源的顶部埋深,然后可计算出均匀异常体的密度值,此反演方法计算快捷,计算量小,可广泛应用于位场数据处理及解释中。 3.本文利用加入了3%高斯噪声的余弦球冠模型数据,反演计算异常源顶部埋深。根据计算的反演数据与理论数据的拟合误差及预测深度与理论模型深度的拟合误差可见,此方法可有效的反演余弦球冠模型的顶部埋深。文中同时给出了不满足条件的反演参数进行反演计算,结果表明,当参数设定范围在视密度及深度不满足线性关系时,反演结果严重失真,观测与预测数据拟合误差及反演与理论模型拟合误差都显著增加。
【学位授予单位】:吉林大学
【学位级别】:博士
【学位授予年份】:2013
【分类号】:P631

手机知网App
【相似文献】
中国期刊全文数据库 前10条
1 朱介寿;骆循;;地球内部三维地震反演方法及其进展[J];国际地震动态;1987年02期
2 单式灶;冯英俊;范祯祥;;二维波动方程速度反演的优化方法[J];石油地球物理勘探;1991年06期
3 黄花香,吴战培,熊艳,王佳川,冯凯;速度反演技术在四川盆地储层预测中的应用[J];石油物探;2002年04期
4 姜彦南,阮帅,王绪本,高永才;核磁共振地下水勘查的模拟退火反演[J];物探化探计算技术;2005年02期
5 方颖,江在森,牛安福;川滇菱形块体东边界地壳形变研究[J];大地测量与地球动力学;2005年03期
6 周仕勇;陈晓非;;近震源破裂过程反演研究——Ⅱ.9.21中国台湾集集地震破裂过程的近场反演[J];中国科学(D辑:地球科学);2006年01期
7 王云专;杨立伟;李素华;;剩余时差校正及泊松比反演[J];地球物理学进展;2006年01期
8 范冬娟;张韶华;;高光谱影像反射率反演方法的研究[J];海洋测绘;2006年03期
9 胡泽勇;程国栋;谷良雷;李茂善;马耀明;;青藏铁路路基表面太阳总辐射和温度反演方法[J];地球科学进展;2006年12期
10 王艳秋;刘旭春;潘华盛;;基于GPS技术的大气可降水反演方法[J];气象科技;2007年03期
中国重要会议论文全文数据库 前10条
1 游然;卢乃锰;李小青;;风云三号微波降水反演[A];第二届微波遥感技术研讨会摘要全集[C];2006年
2 龙桂华;李小凡;张美根;;波动方程多尺度拟牛顿信赖域反演方法的研究[A];中国地球物理学会第二十三届年会论文集[C];2007年
3 徐丽;胡雄;程永强;郭商勇;闫召爱;Xinzhao Chu;;钠测风测温激光雷达探测原理及大气参数反演方法研究[A];第十二届全国日地空间物理学术讨论会论文摘要集[C];2007年
4 闫永利;马晓冰;底青云;孔祥儒;高立兵;陈刚;;层状介质二维电阻率扰动反演方法[A];中国科学院地质与地球物理研究所二○○四学术论文汇编·第一卷(地球深部)[C];2004年
5 刘亚亚;毛节泰;刘钧;李峰;;地基微波辐射计遥感大气廓线的BP神经网络反演方法研究[A];第28届中国气象学会年会——S9大气物理学与大气环境[C];2011年
6 李少鹏;印兴耀;张繁昌;印海燕;;AVO地震参数反演方法研究[A];中国地球物理学会第二十四届年会论文集[C];2008年
7 张文飞;;粘弹性层状介质中的动态多参数反演方法[A];面向国民经济可持续发展战略的岩石力学与岩石工程——中国岩石力学与工程学会第五次学术大会论文集[C];1998年
8 眭素文;于长春;;起伏地形条件下无负值视磁化强度反演方法[A];中国地球物理.2003——中国地球物理学会第十九届年会论文集[C];2003年
9 丁继才;常旭;;基于声波方程的井间地震数据WTW反演新思路[A];中国地球物理学会第22届年会论文集[C];2006年
10 刘展;赵文举;相鹏;;多层密度界面的拟BP神经网络反演方法[A];中国地球物理学会第二十三届年会论文集[C];2007年
中国重要报纸全文数据库 前10条
1 郭志武;风云翘楚看巾帼[N];中国气象报;2002年
2 王五平;声波CT为混凝土大坝做体检[N];中国水利报;2005年
3 许崇宝 徐波 宋嗣麒;山东局物测队煤田岩性地震勘探取得突破[N];中煤地质报;2007年
4 记者 王巧然通讯员 李芬;女博士石玉梅科研成果获SEG高度评价[N];中国石油报;2008年
5 师啸;吐哈油田依靠技术进步加快勘探步伐[N];中国石油报;2008年
6 宁津生 陈军 晁定波;获取地球内部信息的间接方法[N];中国测绘报;2002年
7 记者 魏景云;气象卫星透视地下水[N];中国气象报;2001年
8 朱彦彬;非地震专家何展翔榜上有名[N];中国石油报;2003年
9 本报实习记者 郑菲;争分夺秒 天道酬勤[N];中国气象报;2007年
10 彭冬梅;初至波层析反演静校正技术[N];中国石油报;2004年
中国博士学位论文全文数据库 前10条
1 刘银萍;基于Extrapolation Tikhonov正则化算法的重力数据及梯度多分量数据的3D反演方法研究[D];吉林大学;2013年
2 林雪慧;基于实验的杂交反演方法及其在复合材料性能研究中的应用[D];天津大学;2004年
3 谭代明;隧道超前探水全空间瞬变电磁理论及其应用研究[D];西南交通大学;2009年
4 轩义华;倾斜横向各向同性介质(TTI)参数反演方法研究[D];吉林大学;2007年
5 宋宗平;火山岩储层地震预测技术研究[D];成都理工大学;2008年
6 黄靓;混凝土超声波层析成像的理论方法和试验研究[D];湖南大学;2008年
7 唐红;光全散射法颗粒粒径分布反演算法的研究[D];哈尔滨工业大学;2008年
8 韩利芬;基于神经网络的薄板冲压成形中的反演问题研究[D];湖南大学;2006年
9 吴小成;电离层无线电掩星技术研究[D];中国科学院研究生院(空间科学与应用研究中心);2008年
10 张开鹏;结构损伤识别方法研究[D];武汉理工大学;2004年
中国硕士学位论文全文数据库 前10条
1 成琥;基于BP-GA混合算法的波阻抗反演研究[D];成都理工大学;2007年
2 唐小兵;磁异常处理解释系统的设计与开发[D];成都理工大学;2008年
3 王阅兵;利用Yabuki&Matsu'ura反演方法研究日本东北地震和汶川地震的同震滑动模型[D];中国地震局地震预测研究所;2012年
4 李林峡;高精度谱反演方法研究与应用[D];成都理工大学;2012年
5 葛成;核磁共振T_2谱反演方法研究[D];大庆石油学院;2010年
6 王康宁;高分辨率叠后反演方法研究及在塔河某区的应用[D];中国地质大学(北京);2010年
7 周运波;微震监测反演方法研究[D];长江大学;2012年
8 王勇;基于遗传算法的CSAMT反演研究与应用[D];中国地质大学(北京);2006年
9 仇成强;基于SAR图像的土壤含水量反演方法研究[D];电子科技大学;2012年
10 陈科;基于模拟退火的谱反演方法研究[D];中国石油大学;2010年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026